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ABSTRACT 


This  thesis  investigates  the  simulation  perfonnance  of  Multiple-Input  Multiple-Output 
(MIMO)  communication  systems.  Uncoded  narrowband  and  wideband  MIMO  systems 
are  presented  and  simulated.  The  system  is  further  refined  with  the  addition  of  Space- 
Time  Codes  (STC)  and  Channel  State  Information  at  the  Transmitter  (CSIT).  A 
disadvantaged  receiver  is  subsequently  introduced  to  the  system  lacking  the  optimization 
enjoyed  by  the  native  receiver.  Simulation  and  analysis  was  conducted  with  multiple 
modulation  schemes  and  antenna  configurations.  Rayleigh  and  Rician  fading  models  are 
developed  and  simulated  as  the  wireless  channel.  Receiver  performance  results  based  on 
MATLAB  simulation  are  compared  and  presented. 
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EXECUTIVE  SUMMARY 


Multiple -Input  Multiple-Output  (MIMO)  communications  systems  are  fundamentally 
different  from  earlier  wireless  systems.  In  MIMO  systems,  multiple  antennas  are  used 
for  transmission  and  reception,  thus  multiple  communication  paths  or  channels  are 
developed.  These  channels  can  then  be  utilized  to  increase  the  data  rate  or  reliability  of 
the  system.  Data  rate  is  increased  by  sending  different  symbols  over  multiple  channels  at 
the  same  time.  On  the  other  hand,  reliability  is  increased  by  sending  multiple  copies  of 
the  same  symbol,  or  similar  redundancy,  over  the  multiple  channels. 

This  ability  of  MIMO  systems  to  provide  increased  reliability  and  data  rate 
without  requiring  an  increase  in  bandwidth  or  power  has  resulted  in  MIMO  schemes 
being  adopted  in  standards  such  as  IEEE  802.1  In,  IEEE  802.16,  and  Long  Term 
Evolution  (LTE)  systems,  representing  most  of  the  WLAN,  WMAN,  and  4G  mobile 
telephone  industry  [1],  Our  adversaries  use  conventional  wireless  devices,  and  we  can 
expect  them  to  use  these  MIMO  devices  also,  as  they  grow  in  market  dominance.  We 
benefit  from  a  thorough  understanding  of  our  adversaries’  tools.  Furthermore,  MIMO 
systems  have  the  potential  to  greatly  improve  signal  reception  and  reliability  for  friendly 
forces.  As  the  military  continues  to  develop  and  further  rely  on  high  data  rate,  highly 
mobile  communications,  signal  reception  and  reliability  is  often  found  to  be  the  weakest 
link.  Mobile  networking,  in  particular,  relies  heavily  on  reliable  delivery,  otherwise 
routing  protocols  degenerate  and  retransmissions  swamp  the  network.  By  developing  and 
utilizing  a  robust  MIMO  system,  decision  makers  can  be  assured  that  timely  and  critical 
communications  will  get  through.  This  work  seeks  to  expand  that  understanding. 

Specifically,  this  thesis  evaluates  the  bit  error  rate  performance  of  various  MIMO 
system  configurations  into  which  a  second,  disadvantaged  receiver  is  introduced.  The 
analysis  begins  with  a  basic  two  transmit  and  two  receive  antenna  MIMO  system 
utilizing  maximum  likelihood  decoding.  The  system  is  simulated  in  MATLAB 
employing  BPSK,  QPSK,  16QAM  and  64QAM  over  both  Rician  and  Rayleigh  channel 
models  in  order  to  demonstrate  and  analyze  the  bit  error  rate  performance.  The  system  is 
improved  with  the  application  of  the  Alamouti  space-time  coding  scheme  and  simulated. 


xv 


A  comparison  of  the  uncoded  and  coded  systems  is  conducted.  Specifically,  when 
comparing  a  2x2  MIMO  system  utilizing  QPSK  to  that  of  Alamouti  space  time  coded  2x2 
MIMO  system  utilizing  QPSK,  we  see  that  the  coded  system  provides  the  same 
performance  while  requiring  10  dB  less  power.  The  performance  improvement  seen  is 
due  to  the  diversity  gain  provided  by  the  Alamouti  scheme.  Additionally,  the 
orthogonality  of  the  coding  sequence  transforms  the  decision  at  the  receiver  from  a  vector 
operation  into  a  scalar  operation,  reducing  the  receiver  complexity. 

Next,  the  system  is  further  refined  with  the  addition  of  channel  state  information 
at  the  transmitter.  This  infonnation  allows  the  transmitter  to  optimize  its  transmission 
according  to  the  channel  conditions,  further  improving  bit  error  performance.  Results 
show  that  a  system  with  channel  state  information  needs  12  dB  less  power  than  the 
equivalent  uncoded  system  for  equal  perfonnance.  In  addition,  knowledge  of  the  channel 
provides  the  receiver  with  an  array  gain  equal  to  2  dB,  resulting  in  increased  performance 
over  that  of  the  space-time  coded  system  without  channel  state  infonnation. 

Into  this  robust  system,  a  second  receiver  is  introduced  that  does  not  enjoy  the 
same  optimization  as  the  intended  receiver.  This  would  be  important  whenever  there  are 
multiple  receivers  attempting  to  receive  the  same  transmission  at  distinct  locations,  as 
might  be  the  desired  configuration  in  a  wireless  network  employing  broadcast 
transmissions.  The  performance  was  evaluated  over  both  Rayleigh  and  Rician  fading 
channels  with  various  antenna  configurations.  In  the  case  of  a  Rayleigh  channel,  when 
the  number  of  receive  antennas  is  increased  to  eight,  the  perfonnance  is  3  dB  worse  than 
the  intended  receiver.  In  a  Rician  channel,  where  the  line  of  sight  path  has  equal  power 
to  that  of  the  scattered  paths,  results  show  that  when  the  receiver  employs  four  antennas 
the  performance  increases  to  2  dB  better  than  the  intended  receiver.  Finally,  when  the 
disadvantaged  receiver  is  simulated  at  a  location  receiving  more  signal  power  via  line  of 
sight  from  the  transmitter  with  two  receive  antennas,  the  second  receiver  perfonns  as  well 
as  the  intended  receiver.  Thus,  depending  on  the  configuration  of  the  disadvantaged 
receiver  it  can  achieve  an  equal  or  superior  bit  error  rate  performance  than  the  intended 
receiver. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Multiple-input  Multiple-output  (MIMO)  communications  systems  are 
fundamentally  different  from  earlier  wireless  systems.  As  their  name  implies,  MIMO 
systems  are  designed  to  have  additional  antennas  at  the  transmitter  and  receiver, 
providing  spatial  diversity  not  available  to  single  antenna  systems.  Thus,  information  can 
be  sent  and  received  over  multiple  channels. 

This  gives  the  communications  engineer  greater  opportunity  for  robust  system 
designs.  One  such  example  is  the  capability  of  a  MIMO  system  to  transmit  at  a  higher 
rate  and  more  reliably  by  using  the  stronger  (less  faded)  channels  more  and  the  weaker 
(more  faded)  channels  less.  As  it  is  discussed  in  this  thesis,  that  implies  the  transmission 
is  optimized  for  the  location  of  one  receiver,  which  is  unfortunate  in  the  case  of  broadcast 
transmissions.  This  research  seeks  to  determine  how  a  second  receiver  would  perform 
given  that  the  transmission  is  optimized  for  the  location  of  the  first  receiver. 

B.  OBJECTIVES  AND  METHODOLOGY 

The  goal  of  this  work  was  to  investigate  the  performance  of  a  second  receiver 
when  presented  into  an  existing  MIMO  system  optimized  for  the  first  receiver.  The 
research  approach  takes  a  basic  system,  integrates  techniques  to  improve  performance 
and  concludes  with  the  introduction  of  a  second  receiver  into  the  system. 

Development  begins  with  the  simulation  of  a  simple  uncoded  narrowband  MIMO 
system  over  Rician  and  Rayleigh  fading  channels.  To  which  the  Alamouti  space-time 
coding  scheme  is  applied  and  simulated.  Then,  the  application  of  Orthogonal  Frequency 
Division  Multiplexing  (OFDM)  is  investigated.  The  system  is  completed  by  introducing 
Channel  State  Infonnation  at  the  Transmitter  (CSIT),  where  the  result  is  a  more  robust 
communications  link.  Finally,  a  second  receiver  is  introduced  to  the  optimized  system 
and  the  perfonnance  is  simulated.  Various  scenarios  are  then  investigated  in  an  effort  to 
improve  the  disadvantaged  receiver  performance. 
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C.  BENEFITS  OF  STUDY 

MIMO  systems  are  becoming  increasingly  commonplace,  including  IEEE 
802.1  In,  IEEE  802.16,  and  Long  Tenn  Evolution  (LTE)  systems,  which  represent  most 
of  the  WLAN,  WMAN,  and  4G  mobile  telephone  industry  [1].  Our  adversaries  use 
conventional  wireless  devices,  and  we  can  expect  them  to  use  these  MIMO  devices  also, 
as  they  grow  in  market  dominance.  We  benefit  from  a  thorough  understanding  of  our 
adversaries’  tools.  This  work  seeks  to  expand  that  understanding.  Furthermore,  MIMO 
systems  have  the  potential  to  greatly  improve  signal  reception  and  reliability  for  friendly 
forces.  As  the  military  continues  to  develop  and  further  rely  on  high  data  rate,  highly 
mobile  communications,  signal  reception  and  reliability  is  often  found  to  be  the  weakest 
link.  Mobile  networking,  in  particular,  relies  heavily  on  reliable  delivery,  otherwise 
routing  protocols  degenerate  and  retransmissions  swamp  the  network.  By  developing  and 
utilizing  a  robust  MIMO  system,  it  is  highly  likely  that  timely  and  critical 
communications  will  get  through. 

D.  ORGANIZATION 

The  thesis  is  organized  into  five  chapters.  Chapter  II  introduces  the  MIMO 
system.  It  provides  a  brief  background  into  methods  used  throughout  the  thesis.  These 
include  channel  models,  OFDM-MIMO,  maximum  likelihood  (ML)  detection,  space  time 
coding  and  channel  state  information.  Chapter  III  discusses  the  analysis  and  simulation 
of  MIMO  receivers  that  involve  various  configurations  of  the  methods  introduced  in  the 
previous  chapter.  Chapter  IV  examines  each  of  these  configurations  and  provides  a 
comparison  of  the  results.  Chapter  V  provides  a  summary  and  gives  recommendations 
for  future  study. 
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II.  OVERVIEW  OF  MIMO  COMMUNICATIONS 


This  chapter  provides  background  on  the  techniques  used  in  designing  and 
analyzing  MIMO  systems.  Specifically,  this  chapter  introduces  a  lowpass  equivalent 
MIMO  system  model  upon  which  the  remaining  topics  will  build.  Rayleigh  and  Rician 
fading  channels  are  then  discussed,  providing  a  framework  for  future  analysis.  Next, 
OFDM  implementation  utilizing  the  Discrete  Fourier  Transform  (DFT)  is  presented. 
Maximum  Likelihood  (ML)  detection  at  the  receiver  is  introduced.  Then,  space-time 
codes  are  developed  using  the  Alamouti  scheme.  Finally,  channel  state  information  is 
defined  for  application  at  the  transmitter. 

A.  LOWPASS  EQUIVALENCY 

Wireless  channels  often  have  different  spectral  characteristics  than  that  of  the 
information  signal.  Therefore,  before  entering  the  channel  it  is  often  necessary  to  match 
the  spectral  characteristics  of  the  signal  to  that  of  the  channel.  In  other  words,  the 
information  signal  is  at  baseband,  a  lower  frequency,  and  needs  to  be  converted  to 
passband,  a  higher  frequency  for  transmission.  At  the  receiver,  the  reverse  takes  place, 
converting  from  passband  to  baseband.  Thus,  when  analyzing  or  designing 
communication  systems  it  is  common  practice  to  represent  the  system  in  lowpass 
equivalent  form.  Representing  the  system  in  this  manner  greatly  simplifies  analysis  since 
the  dependency  on  the  carrier  frequency  is  eliminated.  Additionally,  signal  processing 
limitations  require  that  we  manipulate  the  information  signal  at  lower  frequencies  in 
order  to  maintain  reasonable  sampling  rates.  This  section  will  develop  an  equivalent 
lowpass  fonn  for  bandpass  signals  [2],  [3]. 

Given  a  bandpass  signal  g  (7),  whose  Fourier  transform  is  G  ( f ) ,  then  the 
analytic  signal  is  defined  as  [2] 

g+{t)  =  g(t)  +  M{t)  (2-1) 
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where  g+  (t)  denotes  only  positive  frequency  components  are  present  and  g{t) 

represents  the  Hilbert  Transform  of  g(i) ,  i.e.,  g(/)  =  —  *g(/),  where  the  asterisk 

Kt 

denotes  convolution. 

Taking  the  Fourier  transform  of  Equation  (2.1)  yields  [1] 

G.  (/)  =  G(/) + j(-js  gn  (/)  G(/))  (2.2) 

Equation  (2.2)  can  be  simplified  by  using  the  property,  sgn(/)  =  2w(/)-l,  where 
u  (/)  is  the  unit  step  function,  to  get 

G+(f)  =  2G(f)u(f)  (2.3) 

From  Equation  (2.3)  it  is  clear  that  the  analytic  signal  represents  the  bandpass  signal  at 
twice  the  amplitude  and  consisting  of  positive  frequencies  only.  Next,  let  the  lowpass 
equivalent  of  g(t) ,  denoted  as  g,  (t) ,  be  defined  as  [2] 

gl(t)=g+(t)eJ2*"  (2.4) 

Taking  the  Fourier  transform  and  substituting  Equation  (2.3)  results  in 

G,(f)  =  2G(f  +  f')u(f+fc)  (2.5) 

Equation  (2.5)  shows  in  the  frequency  domain,  the  lowpass  equivalent  signal  or 
complex  envelope  is  just  a  scaled  and  shifted  version  of  the  components  of  the  bandpass 
signal  that  have  positive  frequencies.  The  positive  frequency  components  contain  all  the 
information  necessary  to  reconstruct  the  original  signal.  In  this  case,  it  is  scaled  by  two 
and  shifted  from  the  carrier  frequency  to  the  zero  frequency  or  baseband.  This  result  is 
illustrated  in  Figure  1. 
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Figure  1.  A  bandpass  signal  (a)  and  the  lowpass  equivalent  (b)  (After  [2]). 


In  the  time  domain 

&/(0  =  (&(  t)  +  Ji{t))e~j2’cfc‘  (2-6) 

Rearranging  and  isolating  the  real  part  of  Equation  (2.6)  produces  a  representation  of  a 
bandpass  signal  in  terms  of  its  lowpass  equivalent  [2] 

g(t)=Re[g,(t)e72^]  (2.7) 

Additionally,  the  lowpass  equivalent  signal  can  be  represented  by  the  in-phase 
and  quadrature  components  of  the  baseband  signal  as  [2] 

s,(t)=gi(t)+jgq{t) 

Substituting  Equation  (2.8)  into  Equation  (2.6)  yields  [2] 

g,  (0  =g(t)cos(2nfct)+i(t)sm(2nfct)  and 
gq(t)  =  ^(t)cos(2xfct)- g(t)sm(2xfct). 

Solving  for  g(7)  and  yields  [2] 

g  (0  =  g,  {t)cos(2nfct)  -  gq  ( t )  sin  (2 nfj) 

3(0  =  g9  (0 cos (2 jufct)  +  g,  (t)&in(2xfct) 

or,  in  polar  fonn,  [2] 


(2.8) 


(2.9) 


(2.10) 
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S{t)  =  rg(t)cOs(27Tfct  +  0g(  f)) 


(2.11) 


where  ra  (t)  and  6a  (t  )  represent  the  envelope  and  phase  of  g(t) . 

Thus,  we  observe  that  Equation  (2.10)  defines  the  modulation  process  while 
Equation  (2.9)  defines  the  demodulation  process  and  original  signal  can  be  completely 
described  in  lowpass  equivalent  fonn  given  by  Equation  (2.7).  This  conclusion  is  shown 
in  Figure  2.  [2], 


(*) 

Figure  2.  A  modulator  (a)  and  demodulator  (b)  (After  [2]). 

Next,  we  develop  how  the  energy  of  the  signal  is  represented  when  the  signal  is 
expressed  in  lowpass  equivalent  form. 
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(2.12) 


Given  that  the  energy  of  a  signal  g(t)  is  defined  as  [2] 


Es  =  Jk(0f* 


applying  Parseval's  Theorem  yields  [2] 

oo  oo 

ES  =  JkWf*=  JIG(/)|V- 

-00  -00 

Since  G(f)  =  G+(f)  +  G_(f)  and  G+(/)G_(/)  =  0  then  [2] 

oo 

E«  =  |lG,(/)+G(/)|V 

—00 

00  00 

=  }|g.(/)|V+||g(/)|V 

—00  —00 

00 

=  2||G,(/)|V 

—00 


(2.13) 


(2.14) 


This  result  confirms  Equation  (2.5):  the  energy  of  the  lowpass  equivalent  signal  is 
double  that  in  the  bandpass  signal.  In  the  remaining  analysis,  signals  and  systems  will  be 
represented  by  their  lowpass  equivalent  form. 


B.  THE  MIMO  SYSTEM  MODEL 


A  MIMO  system  comprised  of  NT  spatially  separated  transmit  antennas  and  NR 
spatially  separated  receive  antennas  is  depicted  in  Figure  3. 

At  the  transmitter,  the  symbol  transmitted  on  antenna  n  is  denoted  sn(t) ,  where 

n  =  l,2,...,NT  and  t  denotes  the  time  variable.  The  equivalent  lowpass  channel  impulse 
response  between  the  receive  antenna  m ,  and  the  transmit  antenna  n ,  is  denoted  as 
hmn  (x;t) .  Where  T  denotes  the  delay  variable. 
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Figure  3.  MIMO  System  Block  Diagram  (After  [2], [3]). 


The  composite  of  the  complex  channel  gains  is  then  given  by  the  NR  x  Nr  matrix, 


H(x;t)  which  is  also  known  as  the  channel  matrix  [2] 


H(t  ;  t) 


/?n(x;t) 

*12  fat) 

*uvrfat) 

*21  fat) 

*22(x;t) 

Knt  fat) 

Kr  2  fat) 

*wj2fat)  ••• 

* NrNt  (x>  t) 

(2-15) 


The  channel  matrix  is  often  described  with  statistical  models  that  accurately 
predict  channel  characteristics.  This  thesis  will  focus  on  Rayleigh  and  Rician  fading 
models,  discussed  in  the  next  chapter. 


The  signal  at  the  mth  receive  antenna,  ym,  is  the  sum  of  the  transmitted  signals 
corrupted  by  the  channel  and  noise,  resulting  in  [2] 


Nt  co 


y, 


,(0  =  Xj^„(x;tK(t-x)Jx 


n— 1  —oo 
NT 


(2.16) 


m  =  1,2...,  Nr 
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where  the  asterisk  implies  convolution.  The  system  can  therefore  be  expressed  as  a 
system  of  linear  equations  as 


Ti(0 

hn(x;t) 

^12  (Ut) 

V(T;t) 

Sl(x) 

y2(  0 

= 

Mu  0 

h22  (u  t) 

k2Nr  fat) 

* 

u(x) 

i 

5? 

.  ^ 

_ 1 

KR2^)  ••• 

s-  • 

1 _ 

(T)_ 

or  more  conveniently  in  matrix  form  as  [2] 

y(0  =  H(x;t)*s(x)  (2.18) 

where  y  represents  the  NR  x  1  received  matrix  and  s  represents  the  Nr  x  1  sent  matrix 

[2]. 


The  result  above  provides  a  baseline  MIMO  channel  model.  To  refine  this  model 
we  make  the  following  assumptions,  taken  largely  from  [2],  that  will  be  used  throughout 
this  thesis: 

1.  Thennal  noise  is  present  at  each  receive  antenna  and  is  modeled  as  additive 
white  Gaussian  noise  (AWGN),  denoted  as  Tf,,.  Further,  the  AWGN  can  be 
characterized  as  independent  identically  distributed  (IID),  each  having  zero 
mean  and  two  sided  power  spectral  density  2 Na . 

2.  The  receiver  utilizes  a  matched  filter,  whose  output  is  sampled  at  the  end  of 
each  symbol  interval. 

3.  The  intended  receiver  has  perfect  channel  estimation.  Though  this  is  not 
practical,  it  can  be  very  accurately  approximated  through  the  use  of  training 
symbols. 

4.  Antenna  spacing  is  sufficient  so  that  the  fading  on  the  sub-channels  are 
uncorrelated. 

Using  these  assumptions,  the  baseline  MIMO  model  can  be  updated  as  shown  in  Figure 
4. 
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Figure  4.  Updated  MIMO  System  Block  Diagram  (After  [2], [3]). 


Therefore,  hmn  (t)  can  be  simplified  to  hmn ,  and  the  composite  of  the  channel  coefficients 
is  given  by  the  N R  x  Nr  matrix,  H  [2] 


H 


‘21 


‘22 


h 

h 


1  Nt 
2Nt 


h 


NRl 


h 


nrnt 


(2.19) 


The  Nrx  1  noise  matrix,  r|  is  given  as  [2] 


Di 

n2 


(2.20) 


The  received  signal,  ym ,  at  each  antenna  is  the  sum  of  the  transmitted  signals  corrupted 

by  the  channel  and  noise,  resulting  in  [2] 

nt 

ym=T,Snhmn+^  TH  =\,2.-,N R  . 

n= 1 

This  can  be  expressed  as  a  system  of  linear  equations 
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(2.21) 


Ti 

hn 

h\i 

1 - 

’  fix  ’ 

y2 

= 

h2l 

h-22 

h 

n2NT 

V 

+ 

fi2 

(2.22) 

i 

KRl 

h 

nNR2 

h 

hnrnt  J 

snt 

or,  more  conveniently,  in  matrix  fonn  as  [2] 

y  =  Hs  +  r| . 

This  result  will  be  used  throughout  the  thesis. 

C.  MULTIPATH  CHANNEL  MODELS 


(2.23) 


Over  a  wireless  channel,  a  radio  frequency  (RF)  signal  can  take  numerous  paths  to 
its  destination.  This  is  due  to  objects  in  the  propagation  path  reflecting,  diffracting,  and 
scattering  the  radio  frequency  signal  as  it  travels  from  the  transmit  antenna  to  the  receive 
antenna.  These  paths  result  in  multiple  copies  of  the  transmitted  signal  arriving  at  the 
receiver  that  vary  randomly  in  amplitude,  phase  and  time.  Collectively,  this  is  called 
multipath  [5],  [7]. 

Variations  in  amplitude  can  produce  signals  that  are  too  weak  to  be  detected. 
Differences  in  phase  produce  signals  that  constructively  and  destructively  interfere. 
Additionally,  deviation  in  arrival  time  often  results  in  intersymbol  interference.  That  is, 
symbols  which  overlap  each  other  in  time.  In  the  frequency  domain,  relative  motion  of 
the  transmitter,  receiver,  or  the  objects  between  them  results  in  a  Doppler  shift 
proportional  to  the  relative  velocity  of  the  objects.  The  significance  of  these  effects  are  a 
function  of  the  transmitted  signal  and  either  the  coherence  time  ( Tc  )  or  bandwidth  (5C). 

That  is,  how  quickly  the  channel  changes  in  time  and  frequency  respectively.  This  allows 
the  effects  of  the  fading  channel  to  be  classified  by  four  types:  slow  fading,  fast  fading, 
flat  fading,  and  frequency  selective  fading.  This  fading  classification,  relative  to  the 
signal  bandwidth  ( Bs )  and  symbol  time  ( Ts ),  is  depicted  in  Table  1.  [2],  [5],  [7]. 

The  coherence  time  of  the  channel,  Tc  depends  on  the  Doppler  spread,  Bn 
according  to  the  relationship  [6] 
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1 


(2.24) 


Tc  — 

4Bd 

where  the  Doppler  spread  is  a  measure  of  the  spectral  expansion  due  to  the  relative 
motion  between  the  transmitter  and  receiver.  If  the  coherence  time  is  short  relative  to  the 
symbol  duration,  then  the  channel  impulse  response  changes  rapidly  over  the  symbol 
period,  causing  distortion.  The  channel  is  said  to  be  fast  fading.  On  the  other  hand,  if  the 
coherence  time  is  large  relative  to  the  symbol  period  then  the  channel  can  be  considered 
static  over  multiple  symbols.  The  channel  is  said  to  be  slow  fading.  Thus,  classification 
of  a  channel  as  either  fast  or  slow  fading  is  detennined  by  how  quickly  the  channel  is 
changing  in  time  and  the  baseband  signaling  utilized  [5]-[7]. 

The  coherence  bandwidth  of  the  channel,  Bc  depends  on  the  multipath  delay 
spread,  TD  by  the  relationship  [6] 

Bc=^r  (2.25) 

‘  D 

where  the  multipath  delay  spread  is  a  measure  of  the  difference  of  arrival  times  between 
the  longest  and  shortest  propagation  paths.  Only  paths  with  significant  energy  are 
considered.  If  propagation  delay  is  large  (i.e.,  large  TD)  with  respect  to  the  symbol 

duration,  then  multiple  symbols  overlap  each  other  at  the  receiver,  causing  distortion. 
The  channel  is  said  to  be  frequency  selective.  Conversely,  if  the  propagation  delay  is 
small  compared  to  the  symbol  duration,  then  each  path  that  the  symbol  takes  arrives  at 
the  receiver  before  the  next  symbol.  The  channel  is  said  to  be  frequency  non-selective  or 
flat  faded.  Classification  of  the  channel  as  either  frequency  selective  or  flat  faded  is 
determined  by  how  quickly  the  channel  is  changing  in  frequency  and  the  baseband 
signaling  utilized  [5]-[7]. 
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Table  1.  Characteristics  of  channel  fading  (After  [6], [7]). 


Fast  Fading 

^  <^D 

T3>% 

Slow  Fa  ding 

Bf  >  Bn 

TS<TC 

Frequency  Selective  Fading 

BS>BC 

T„>Tt 

Frequency  Nou-Selective  Fading  (Flat) 

A  A 

Therefore,  a  channel  can  be  characterized  and  is  often  referred  to  by  these 
metrics.  This  thesis  will  focus  on  two  statistical  models  that  can  be  used  to  represent  a 
frequency  non-selective  multipath  channel:  Rayleigh  fading  and  Rician  fading. 
Additionally,  the  channel  is  considered  slowly  fading  and  the  complex  channel  gain  is 
constant  over  consecutive  symbols. 

1.  Rayleigh  Fading  Model 

When  the  communicating  antennas  do  not  have  a  direct  line  of  sight  (NLOS)  to 
each  other,  the  received  envelope  can  be  modeled  as  a  Rayleigh  random  process  and  is 
referred  to  as  Rayleigh  fading.  In  order  to  model  Rayleigh  fading,  we  utilize  the  sum  of 
two  IID  Gaussian  random  variables  with  zero  mean. 

Let  ZNLOS  represent  the  sum  of  two  IID  Gaussian  random  variables  with  zero 
mean  and  equal  variance 

^nlos  ~  X + jY  (2.26) 

where  X  ~  TV  (0,  a2 ),  Y  ~  N  (0,  a1 )  and  the  magnitude  R  is  given  by 

Rn,os  =  VA2  +  F2  (2.27) 
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In  order  to  relate  Equation  (2.27)  to  a  distribution,  we  introduce  the  joint  cumulative 
distribution  function  (CDF),  defined  as  [9] 

Fxtx2  {xl,x2)  =  ^r[Xl<xl,X2<x2] 


where  the  CDF  is  related  to  the  joint  probability  density  function  (PDF)  by  [9] 

3  FXlX1  (  X\  >  X2  ) 


fxxxSX  i’x2)  =  - 


dx2dx2 


(2.28) 


(2.29) 


Solving  for  the  CDF  yields  [9] 

X,  x2 

Fx,x2  {x1 , x2 )  =  |  |  fXiXi  (xt,x2) dxtdx2 


—oo  —oo 


thus,  expressing  Equation  (2.27)  in  terms  of  the  CDF  results  in 

Fk  (>')=  ||  fxr  (x,y)dxdy 

\lx2+Y 2  <r 


(2.30) 


(2.31) 


For  r>0.  Since  X  and  Y  are  Gaussian  random  variables,  we  may  further  refine 
Equation  (2.31)  as 


'  i  ' 


vV2/rcr~y  Jx2'+y2 < 


2  2 

-*  -y 


||  elc7~  e2a~  dxdy 


(2.32) 


Simplifying  Equation  (2.32)  and  perfonning  a  change  of  variables  results  in 

(  -,2  A 


\-elF 
v  J 


(r) 


(2.33) 


where  the  unit  step  function  u  is  given  by 


!0)  = 


[  1  r  >  0 

I  0  r  <  0 


(2.34) 


Finally,  given  Equation  (2.33),  solving  for  the  PDF  yields 


/*(r)  =  - 


a 


-e2crw(V) 


(2.35) 
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which  is  the  PDF  of  a  Rayleigh  distribution.  Thus,  the  sum  of  two  IID  complex  Gaussian 
random  variables  with  zero  mean  as  given  by  Equation  (2.26),  accurately  model  Rayleigh 
fading  and  provide  the  desired  channel  model  when  NLOS  components  are  present  at  the 
receiver  [2], 


2.  Rician  Fading  Model 


For  channels  with  a  line  of  sight  (LOS)  path  in  addition  to  NLOS,  the  Rician 
model  is  more  appropriate.  However,  due  to  the  LOS  path,  the  Gaussian  random 
variables  no  longer  can  be  modeled  with  equal  mean. 


The  Rician  model  follows  a  Rician  distribution  given  by  [7] 


(2.36) 


where  A  >  0  and  I0  (□)  represents  the  modified  Bessel  function  of  the  first  kind  and  zero 
order  [7]. 


Alternatively,  this  distribution  can  be  described  in  terms  of  the  A'- factor,  which  is 
defined  as  the  ratio  of  the  energy  in  the  specular  (LOS  path)  to  the  energy  in  the  diffuse 
(NLOS  paths),  given  by  [8] 


K  = 


(2.37) 


The  second  moment  or  average  received  power  is  the  total  power  of  the  specular 
and  diffuse  paths,  given  by 

P  =  A2  +  la2  (2.38) 


Solving  for  A2  and  2 a2  and  substituting  Equation  (2.37)  results  in 


A2 

2  a2 


KPr 
K  + 1 
Pr 

K  + 1 . 


(2.39) 
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Now,  the  X-factor  can  be  used  to  represent  a  complex  random  variable  whose 
magnitude  is  a  Rician  random  variable,  given  in  Equation  (2.36),  by  utilizing  Equation 
(2.39)  as  [8] 


—  - 
'K  + 1 


-ere 


je 


(2.40) 


where  the  first  term  is  from  the  specular  path  and  the  second  term  is  from  the  diffuse 
path.  The  diffuse  path  is  equivalent  to  Rayleigh  fading  described  in  the  previous  section. 

Thus,  by  Equation  (2.40),  Rician  fading  can  be  modeled  by  the  sum  of  a  LOS  path 
with  uniform  phase  and  NLOS  paths  as  characterized  by  Equation  (2.26)  scaled  by  the 
X-factor.  It  is  easily  seen  that  if  K  =  0,  Equation  (2.40)  reduces  to  a  zero  mean  complex 
Gaussian,  which  can  be  used  to  model  the  hmn ,  for  complex  channel  gains.  Additionally, 

as  K  becomes  very  large,  Equation  (2.40)  becomes  detenninistic  and  would  be  used  for 
the  non-faded  case  [2]. 

D.  ORTHOGONAL  FREQUENCY  DIVISION  MULTIPLEXING 


A  primary  reason  for  employing  MIMO  systems  is  to  achieve  higher  data  rates 
when  compared  to  single  antenna  systems.  However,  as  the  data  rate  increases,  the 
symbol  time  decreases  relative  to  the  multipath  delay  spread.  As  discussed  in  the 
previous  section,  this  produces  frequency  selective  fading  which  leads  to  intersymbol 
interference  (ISI)  and  deep  fades,  severely  degrading  performance  [10], [1 1]. 

To  mitigate  ISI,  multicarrier  modulation  schemes  have  been  developed  that  divide 
the  large  bandwidth  signal  into  L  smaller  subchannels,  where  the  subchannels  are 

characterized  by  —  □  B(  .  That  is,  the  bandwidth  of  each  subchannel  is  equal  to  the 
L 

signal  bandwidth  ( Bs  )  divided  by  the  number  of  total  subchannels  (L).  Thus,  when  L  is 

selected  to  be  sufficiently  large,  the  bandwidth  of  each  subchannel  is  less  than  the 
coherence  bandwidth  of  the  channel  and  each  subchannel  experiences  flat  fading 
conditions.  The  subchannels  are  then  transmitted  in  parallel  such  that  the  bandwidth  of 
the  subchannels  sum  to  that  of  the  original  signal  without  the  effects  of  ISI.  When  the 
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subchannels  are  separated  in  frequency  and  manipulated  to  be  orthogonal  to  each  other, 
the  multicarrier  modulation  scheme  is  called  orthogonal  frequency  division  multiplexing 
(OFDM)  [3], [10]. 

Let  sk  (?)  represent  the  carrier  signal  on  the  kth  subcarrier  [2] 

sk(t)  =  cos27zfkt  (2.41) 

where  k  =  0,1,2,...,Z -1  and  fk  is  the  center  frequency  of  the  klh  subchannel.  The 

subchannels  are  orthogonal  and  independent  of  the  phase  between  subchannels,  over  a 
symbol  interval  T ,  when  the  integral  of  their  product  over  that  same  interval  is  zero, 
given  by  [2] 

T 

J cos  (2 nfkt  +  <pk )  cos  [in  f jt  +  )  dt  =  0  (2.42) 

o 

where  (f>k  and  ^  represent  the  phase  of  the  kth  and  jth  subchannel  respectively  and  are 
independent  of  f k  and  /. .  Selecting  the  OFDM  subchannel  frequency  spacing  to  be 

Y  and  maintaining  the  subchannel  center  frequencies  to  be  integer  multiples  of  each  other 

(i.e.,  fk+\-fk=~  k  =  0,2,...,Z-2)  results  in  the  peak  of  each  subchannel 
corresponding  to  the  nulls  of  adjacent  subchannels.  Thus,  a  set  of  L  parallel  symbols 
spaced  apart  produces  L  orthogonal  overlapping  signals,  providing  bandwidth 
efficiency  and  minimizing  the  effects  of  ISI  [2]. 

1.  OFDM-MIMO  Transmission 

A  block  diagram  showing  the  transmission  OFDM-MIMO  system  with  two 
transmit  antennas  ( Nr  =2)  is  given  in  Figure  5.  At  the  input,  data  is  converted  from  a 
binary  stream  into  a  set  of  L  parallel  symbols  divided  between  each  antenna  and 
modulate,  where  the  output  of  the  modulator  is  of  the  form  X  [0]  ,X[A-l]  and 
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the  symbol  rate  is  reduced  by  a  factor  of  —  containing  bk  bits.  Modulation  is  typically 

L 

chosen  to  be  consistent  across  subcarriers.  However,  it  is  possible  to  modulate  each 
subcarrier  with  a  different  scheme. 

When  L  is  selected  to  be  an  integer  multiple  of  two,  then  an  efficient  method 
known  as  the  Inverse  Fast  Fourier  Transform  (IFFT)  can  be  used  to  implement  the 
inverse  discrete  Fourier  transform  (IDFT)  given  by 

i  L- 1  jinlk 

j Al\  =  -rlLX[kY  L  /  =  0,l,...,z-l  (2.43) 

\L  k= o 


As  can  be  seen  from  Equation  (2.43),  the  IDFT  takes  the  components  of  each 
subcarrier  and  produces  equivalent  time  domain  symbols  that  have  orthogonal  properties 
mentioned  previously. 


Figure  5.  OFDM-MIMO  Transmission  Block  Diagram  (After  [2]). 

The  next  step  in  the  transmission  process  is  the  addition  of  a  cyclic  prefix  (CP). 
The  CP  is  created  by  replicating  the  data  contained  in  the  last  part  of  the  symbol  and 
superimposing  this  copy  to  the  front  of  the  symbol  as  depicted  in  Figure  6.  The  total 
symbol  duration  is  given  by 

Ttot  =  Tcp  +  Ts  (2.44) 

This  is  done  to  guard  against  the  inherent  delay  of  a  multipath  channel.  As 
introduced  in  previous  sections,  as  the  symbol  propagates  through  the  wireless  channel,  a 
delay  is  incurred  at  the  receiver.  The  cyclic  prefix  is  selected  to  be  of  sufficient  length  as 
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to  prevent  the  delay  of  one  symbol  from  corrupting  the  following  symbol.  In  other 
words,  the  cyclic  prefix  is  used  as  a  guard  interval  to  keep  the  transmitted  symbols 
independent.  Obviously,  in  order  to  do  this,  the  CP  length  must  be  chosen  greater  than 
the  multipath  delay  spread  of  the  channel.  Another  important  aspect  of  the  CP  and  the 
reason  it  is  appended  to  the  front  of  a  symbol  is  to  protect  the  start  point  of  the  symbol 
which  aids  in  bit  decisions  at  the  demodulator  [10]. 


Ter  Tg 

Figure  6.  CP  of  an  OFDM  Symbol  (After  [10]). 


The  symbols  are  then  converted  from  parallel  to  a  serial  sequence,  passed  through 
the  digital  to  analog  (D/A)  converter,  and  unconverted  to  radio  frequency  (RF)  for 
transmission. 

2.  OFDM  Reception 

A  block  diagram  showing  the  OFDM-MIMO  system  reception  process  with  two 
receive  antennas  (NR  =2)  is  given  in  Figure  7.  After  reception,  the  signal  is  down 

converted  and  sent  through  the  analog  to  digital  (A/D)  converter.  The  signal  is  then 
transformed  from  a  serial  sequence  into  parallel  subchannels,  where  the  CP  is  removed 
and  the  discrete  Fourier  transform  (DFT)  is  implemented  with  the  Fast  Fourier  Transform 
(FFT)  which  is  a  Fourier  pair  with  the  IDFT.  Where  the  DFT  is  given  by 

L- 1  -jljtlk 

X[k]  =  ^x[/]e  1  (2.45) 

1=0 

Finally,  the  subchannels  are  demodulated  and  converted  back  into  the  binary  bit 

stream. 
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Figure  7.  OFDM-MIMO  Reception  Block  Diagram  (After  [2]). 

An  OFDM  system  must  maintain  synchronization  and  orthogonality  in  order  to 
accurately  reproduce  the  original  signal.  If  implemented  correctly,  the  CP  will  maintain 
orthogonality  by  preventing  symbol  overlap  as  illustrated  by  Figure  8.  Synchronization 
can  be  maintained  by  allocating  OFDM  symbols  for  this  purpose.  This  thesis  will 
assume  that  the  receiver  has  perfect  synchronization  and  the  orthogonality  of  the  symbols 
is  maintained. 


OFDM  Symbol  A 

B 

a 

OFDM  Symbol  B 

Delay  Spiead 

Figure  8.  ISI  Representation  at  the  Receiver  (After  [10]). 

E.  MAXIMUM  LIKELIHOOD  DETECTION 


Recall  that  the  received  signal,  ym  at  each  antenna  in  a  MIMO  system  is  the  sum 
of  the  transmitted  signals  corrupted  by  the  channel  and  noise.  This  was  expressed  in 
matrix  form  by  Equation  (2.23)  and  is  repeated  here,  y  =  Hs  +  r\ . 

The  corruption  introduced  by  the  wireless  channel  requires  that  the  receiver  must 
statistically  decide  on  which  symbol  was  sent  based  on  what  was  received.  To  make  this 
decision,  the  receiver  utilizes  knowledge  of  the  channel  and  compares  that  to  the  received 
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signal.  For  the  MIMO  model  presented  in  this  thesis  the  optimal  method  to  find  this 
decision  statistic  is  known  as  Maximum  Likelihood  (ML)  detection.  [2] 


Bayes'  rule  forms  the  basis  of  the  ML  detector  and  is  given  by  [9] 

Pr[A2|  A,]Pr[A,] 


Pr[X1|X2]  = 


Pr[X2] 


(2.46) 


Let  Pr[fifj]  represent  the  a  priori  probability  that  a  particular  symbol  vector  s 
was  sent,  where  s  is  selected  from  a  known  signal  space.  Let  Pr[X2]  denote  the  a  priori 
probability  that  a  particular  NR  x  1  decision  statistic  vector  y  was  received. 
Additionally,  let  PrffiGIA]]  indicate  the  a  posteriori  probability  of  y  given  s  was 
observed  [12]. 

Observing  that  sent  symbol  s  is  selected  from  a  finite  set  of  values  and  the 
received  vector  y  is  stochastic,  Equation  (2.46)  can  be  represented  more  accurately  as 


pr(«|y)-/(yl‘)?(,) 
'  /( y) 


(2.47) 


Therefore,  the  detector  needs  to  select  the  sent  symbol  vector  that  maximizes  the 
probability  of  the  sent  symbol  vector,  given  the  received  vector.  Updating  Equation 
(2.47)  results  in  [12] 

/(y  I  s)pr(s) 


s  =  arg  max  - 


f(y) 


(2.48) 


where  s  denotes  the  estimate  out  of  the  detector.  Assuming  that  the  symbols  are  sent 
with  equal  probability,  Equation  (2.48)  becomes  [2] 

s  =  arg  max  /  (y  |  s)  (2.49) 

Since  the  noise  at  the  receiver  is  assumed  to  be  IID  Gaussian  random  variables 
with  zero  mean,  Equation  (2.49)  can  be  represented  as  the  minimization  of  the  distance  to 
the  stationary  vector  y ,  as  [13] 

s  =  argmin|y-Hs|  (2.50) 

S  1 
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Thus,  the  detector  compares  every  possible  symbol  in  the  signal  space  to  the  received 
symbol  and  selects  the  closest.  Alternatively,  Equation  (2.50)  can  be  expressed  in  tenns 
of  each  received  symbol  [2] 


s  =  arg  min  ^ 


m= 1 


y  —  ^h  s 

s  m  mn  n 


n= 1 


(2.51) 


F.  SPACE  TIME  CODING 

Space-time  coding  implies  multiple  antennas  either  at  the  transmitter,  receiver,  or 
both.  They  exist  to  exploit  transmit  diversity  in  order  to  increase  channel  capacity  and 
reliability.  This  thesis  concentrates  on  the  Alamouti  space-time  code,  introduced  in  [14]. 

Figure  9.  depicts  the  scheme  when  applied  to  a  MIMO  system  consisting  of  two 
transmitters  and  two  receivers.  At  time  t  symbols  s,  and  s2  are  transmitted 

simultaneously.  At  time  t  +  T ,  the  negative  complex  conjugate  of  s2  and  the  complex 
conjugate  of  5,  are  transmitted.  The  result  is  a  sequence  of  signals  shown  in  Table  2. 
[14]. 


Hi <*+!> 


Figure  9.  Alamouti  Space-Time  Coding  Scheme  (After  [14]). 

Assuming  that  the  channel  remains  constant  over  at  least  two  consecutive 
symbols,  the  received  signals  can  then  be  expressed  as  [14] 
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(2.52) 


yl(t)  =  \xsl+hl2s2+\\l{t) 
y2(t)  =  h2lsl+h22s2+r\2(t) 
y1(t+T)  =  -hl  xs*  +  hus*  +  pt  (t + T) 
y2(t+T)  =  -h2ls*  +  h22s*  +  r|  2(t  +  T) 

where  yx(t)  and  y,  (t  +  T)  correspond  to  the  received  signals  at  antenna  one  at  time  t 

and  t  +  T  respectively.  Additionally,  p,  (t)  and  p,  (t  +  T)  represent  the  noise 

components  at  antenna  one  at  time  t  and  t  +  T  respectively.  As  mentioned  in  the 
previous  section,  hmn  indicates  the  channel  gain  to  antenna  m  from  antenna  n  . 

Table  2.  Alamouti  scheme  transmission  sequence  (After  [14]). 


Time  t 

Time  t+T 

Tr.in.smit  Antenna  1 

* 

• 

Transmit  Antenna  2 

• 

* 

The  received  signals  are  combined  and  compared.  Resulting  in  symbol  estimates 
given  by  [14] 


S  =  Kyx  (0 + \2yx  (t+T)  +  h2*y2  (t) + Kji  (t+T) 

K  =  hn*y  (t)-Kyx{t+T)  +  Ki  *y2  (0  -  Kyi  (t+T) 


Substituting  Equation  (2.52)  into  Equation  (2.53),  yields 

R  =  hn  ( Vl  +  ^12^2  +  h,  (0)  +  hl2  (~hnS2  +  hl2Sl  +  R  {t  +  TS  ) 
+/?2i  (h2xsx  +  h22s2  +  p2  (t))  +  h22  ^~h2l  s2  +  h22  sx  +  p2  {t  +  T)  j 

P2=*L2*(Vl  +  V2  +rll  (0)-^l  (~\lS2  +hl2Sl  +rll  (^  +  :r)*) 
+^22  (^2  A  ^22S2  ^2  (0)  _  K  ^21  ^2  ^22  ^1  ^2  j 

simplifying 
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(2.55) 


Q=(\hn\2+K\2  +  \hn  |2+ 1*22 12)^1 

~>rh\  I  T]i  (/)  +  ^12^1l  (/  ^21  ^2  (0  ^22^2  (/ 

s = (hi  r + i*i2  r  +[*21  r +i^22  r)^2 

-^n'Hi  {'  +  T)*  +  hu\  (t)  -  h2lr\2  (t  +  T)*  +  h22r \2  ( t  +  T ) 


From  Equation  (2.55)  it  is  clear  that  the  combined  symbols  are  heavily  affected  by  the 
channel  gains.  Therefore,  the  Alamouti  scheme  provides  transmit  diversity  such  that  the 
effects  of  multipath  fading  are  minimized.  The  symbols  EJ  and  E,  are  sent  to  the 
maximum  likelihood  detector,  discussed  in  the  previous  section,  and  results  in  symbol 
estimates  ,v,  and  s[ . 


Another  benefit  of  the  Alamouti  scheme  is  that  it  achieves  full  diversity  while 
maintaining  low  complexity  at  the  receiver.  This  is  due  to  the  orthogonality  of  the 
sequence  of  signals,  shown  in  Table  2.  The  coding  matrix  is  given  by  [14] 


G  = 


■h 

S2 


(2.56) 


where  G  is  the  Alamouti  generator  matrix  for  a  2x2  MIMO  system.  It  can  be  seen  that 
the  row  vectors  are  orthogonal 

Vv/=0  (2.57) 

where  y,  =  [Es’,  -s2*]  and  v2  =  [,v2  .v,"] .  Due  to  this  orthogonality,  the  maximum 

likelihood  detector  no  longer  needs  to  compare  each  possible  transmitted  symbol  to  the 
one  received.  Thus,  the  complexity  of  the  ML  detector  is  reduced  from  M 2  to  2 M 
possibilities,  where  M  represents  the  number  of  possible  received  symbols  in  a  given 
modulation  scheme  (i.e.,  M-QAM).  It  is  clear  that  this  will  significantly  reduce 
computational  time  requirements  [2]. 


G.  CHANNEL  STATE  INFORMATION  AT  THE  TRANSMITTER 


When  the  transmitter  has  access  to  channel  characteristics  by  means  of  a  feedback 

loop  or  another  method,  it  is  said  to  know  the  channel  state.  With  channel  state 
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information  at  the  transmitter  (CSIT),  capacity  can  be  increased  by  allocating  more 
power  to  antennas  with  favorable  channels  [2]. 

Let  the  channel  gain  matrix  H  be  known  to  the  transmitter  and  receiver.  The 
singular  value  decomposition  (SVD)  of  H  is  given  by 

H  =  UEV"  (2.58) 

where  U  is  an  NRxNR  matrix,  V  is  a  NT  x  Nr  matrix  and  both  are  unitary,  that  is 

\}H  U  =  1^  and  \H\  =  1^  .  I  is  a  NR  xNT  diagonal  matrix  of  the  singular  values  of  the 

channel,  denoted  by  yJZt  ,  where  i  is  detennined  by  the  rank  of  H  denoted  by  r  and 

given  by  r  <min(iVr,jVfl)  .  A“  denotes  the  Hennitian  operator,  AH  =  (A  j  [4]. 

Substituting  Equation  (2.58)  into  the  received  signal  vector  given  by  Equation 
(2.23),  yields 

y  =  USV"s+'n  (2.59) 

Defining  linear  transformations  on  the  received  vector  as 

y  =  U^y  (2.60) 

the  sent  vector  as 

s  =  V"s  (2.61) 

and  noise  vector  as 

ri  =  U"r|  (2.62) 

then  applying  these  linear  transformations  to  Equation  (2.59),  yields 

y  =  U"(UEV"s  +  Ti) 

=  U"UEV//Vs  +  U"'n  (2.63) 

=  Es  +  r| 

where  rj  is  distributed  identically  to  r| .  When  the  transformation  of  Equation  (2.61)  is 
applied  to  the  channel  input  s ,  it  is  often  called  transmit  precoding.  When  the 
transformation  of  Equation  (2.60)  is  applied  to  the  output  y  ,  it  is  referred  to  as  receiver 
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shaping.  From  Equation  (2.63),  it  is  clear  that  the  SVD  operation  allows  the  channel  to 
be  represented  by  rank,  r  parallel  channels,  as  is  illustrated  in  Figure  10.  [2]. 


Figure  10.  Equivalent  MIMO  model  (After  [4]). 


This  thesis  will  focus  on  the  ability  of  the  transmitter  to  reduce  errors  by 
allocating  power  to  the  individual  antennas  based  on  the  dominant  eigenvalues  or 

singular  values  of  the  channel  .  This  process  is  repeated  every  symbol  period  to 

reduce  the  probability  of  errors.  If  the  channel  state  is  unknown,  then  the  power  is 
distributed  evenly  among  all  antennas  [2], [3]. 

Given  the  following  total  power  constraint  [4] 

P^P:  (2.64) 


where  Pf  is  the  power  allocated  to  antenna  i  and  P  is  the  total  power.  Let  the  channel 


capacity  be  given  as  [4] 


nt  f  p2  h 

c=sf>g2 

1=1  V  a  ) 


(2.65) 
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where  C  is  the  channel  capacity,  B  represents  the  bandwidth,  and  a1  is  the  noise  power 
of  each  element  of  rj .  Optimization  of  the  channel  capacity  in  terms  of  the  antenna 
power  allocation  is  given  by  [4] 


P  =  max 


—  ,  0 

V  J 


(2.66) 


where  //  denotes  a  constant  that  is  adjusted  until  Equation  (2.64)  is  satisfied.  If  the 
channel  gain  is  low,  then  reliable  communication  is  unlikely  and  the  transmitter  allocates 
less  power  to  that  antenna.  Iterations  are  conducted  in  this  manner  until  the  transmitter 
deems  the  remaining  channels  adequate,  resulting  in  optimal  power  distribution  [4], [5], 


H.  SUMMARY 


This  chapter  introduced  the  MIMO  system  model  which  forms  the  basis  for 
follow-on  analysis.  The  effects  of  multipath  fading  on  a  wireless  channel  were 
introduced  and  two  statistical  models  defined  to  aid  in  analysis.  The  fundamentals  of 
orthogonal  frequency  division  multiplexing  (OFDM)  were  presented  and  a  MIMO- 
OFDM  system  model  was  introduced  utilizing  the  discrete  Fourier  transform  (DFT).  An 
optimal  detection  scheme  in  the  fonn  of  Maximum  likelihood  (MF)  at  the  receiver  was 
discussed.  Then  the  Alamouti  Space-time  coding  method  for  a  2x2  MIMO  system  was 
introduced.  Finally,  channel  state  information  is  defined  for  application  at  the  transmitter 
where  a  power  allocation  process  was  developed.  MIMO  systems  containing  various 
aspects  presented  in  this  chapter  are  developed  and  simulated  in  the  following  chapter. 
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III.  ANALYSIS  AND  SIMULATION 


This  chapter  provides  bit  error  performance  analysis  and  simulation  of  the  two 
transmit  and  two  receive  antenna  baseband  MIMO  model  presented  in  the  last  chapter. 
The  analysis  begins  with  an  uncoded  narrowband  system  and  progresses  by 
implementing  the  various  schemes  introduced  in  the  previous  chapter,  culminating  with  a 
system  utilizing  CSIT.  To  this,  a  second  system  is  added  and  analyzed  where  this  second 
receiver  enjoys  none  of  the  benefits  intended  for  the  primary  receiver.  Each  system  was 
coded  and  simulated  in  MATLAB.  The  source  code  can  be  requested  from  the  manager 
of  the  Communications  Research  Laboratory  (CRL),  at  the  Naval  Postgraduate  School. 

A.  SIMULATION  METHODOLOGY  AND  ASSUMPTIONS 

The  following  assumptions  are  made  when  conducting  the  simulations,  which  are 
described  in  detail  in  the  previous  chapter. 

Given  a  fixed  wireless  communication  system  comprised  of  two  transmit  and  two 
receive  antennas,  where  each  antenna  is  spatially  and  angularly  independent  from  the 
other.  The  information  transmitted  from  each  antenna  uses  the  same  modulation,  either 
MPSK  or  MQAM,  selected  from  a  signal  space  as  illustrated  in  Figure  11.  The 
probability  of  bit  error  (Ph)  for  each  scheme  is  inversely  related  to  the  distance  between 

closest  symbols  and  is  a  measure  of  the  number  of  bits  in  error  per  bit  sent  over  a  given 
time  interval.  Furthermore,  the  modulation  utilizes  Gray  coding,  in  that  adjacent  symbols 
differ  by  only  one  bit.  Both  signals  are  transmitted  with  equal  power,  split  among  the 
two  antennas.  The  total  received  power  from  all  receive  antennas  per  symbol  time  is 
given  by  [19] 

z£[k.r]=>  <3->) 

n 

Due  to  the  statistical  independence  of  |s„|,  Equation  (3.1)  can  be  expressed  for 
each  individual  value  as 
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(3.2) 


1 

J  '  Nr 

In  the  remaining  analysis  and  simulations,  for  each  modulation  scheme,  Eb  is 
determined  and  held  constant  while  Na  is  adjusted  according  to  the  simulated  Eh/Na  . 
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Figure  11.  Signal  Space  for  (a)  BPSK,  (b)  QPSK,  (c)  16QAM,  (d)  64QAM. 


The  channel  is  slow  frequency  non-selective  fading  and  is  modeled  as  having  an 
AWGN,  Rayleigh,  or  Rician  distribution.  Additionally,  the  channel  gains  hmn  are 
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identically  distributed  and  statistically  independent  from  each  other  and  their  sum  is 
normalized  to  unity,  given  by  [2] 


I#, 


=1 


(3.3) 


Due  to  the  statistical  independence  of  hmn ,  Equation  (3.3)  can  be  expressed  for 
each  individual  value  as 


E 


1 


(3.4) 


Thus,  the  average  total  energy  per  bit  received  is  the  transmitted  energy  per  bit 
per  transmit  antenna.  In  other  words,  increasing  the  number  of  transmit  or  receive 
antennas  reduces  the  channel  gains  of  the  system.  This  may  seem  counter  intuitive,  since 
we  would  expect  that  a  system  employing  a  greater  number  of  antennas  would  result  in 
more  received  power  and  thus  larger  energy  per  bit.  However,  this  assumption  is 
consistent  with  the  literature,  namely  [2],  [3],  [5],  and  will  be  utilized  in  this  thesis. 


Synchronization  is  maintained  perfectly  and  the  channel  is  known  to  the  receiver, 
but  unknown  to  the  transmitter.  The  receiver  noise  is  modeled  as  independent  and 
identically  distributed  complex  white  Gaussian  noise  processes,  each  having  zero  mean 
and  power  spectral  density  of  2Na .  Detection  is  accomplished  by  utilizing  a  maximum- 

likelihood  detector,  which  provides  a  best  case  for  probability  of  error.  With  the 
exception  of  equal  power  transmission  and  channel  state  infonnation  at  the  transmitter, 
which  will  be  addressed  specifically,  these  assumptions  will  be  used  in  the  proceeding 
analysis  and  simulation  [2]. 


B.  UNCODED  NARROWBAND 


The  first  task  is  to  evaluate  the  performance  of  an  uncoded  MIMO  system,  as 
depicted  in  Figure  12.  The  entire  system  is  simulated  in  MATLAB  utilizing  equivalent 
baseband  form  in  the  discrete  time  domain.  Two  million  symbols  are  generated  and 
modulated  with  either  BPSK,  QPSK,  16QAM  or  64QAM.  This  is  simulated  for 
increasing  values  of  Eh/Na  over  the  range  of  30  dB. 
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In  these  simulations,  the  transmitter  employs  two  antennas  with  which  to  transmit 
the  modulated  symbols  with  equal  power.  The  wireless  channel  is  modeled  as  slowly 
fading  where  the  signal  period  is  less  than  the  coherence  time.  The  channel  gains  are 
simulated  according  to  the  AWGN,  Rayleigh  and  Rician  models  as  presented  in  Chapter 
II,  under  the  assumptions  given  in  the  previous  section.  At  the  receiver,  a  maximum 
likelihood  decision  statistic  is  simulated  in  which  the  symbol  with  the  smallest  Euclidian 
distance  is  selected.  This  symbol  estimate  is  then  demodulated  according  the  appropriate 
scheme  and  compared  to  the  transmitted  signal.  Decision  errors  are  then  tallied  and 
calculated  over  each  iteration  of  Eb/Na  resulting  in  the  following  bit  error  plots. 


Figure  12.  MIMO  System  Block  Diagram  (After  [2], [3]). 


In  addition,  for  comparison  and  to  validate  the  accuracy  of  the  program,  an 
AWGN  channel  is  simulated  and  the  ideal  bit  error  probability,  Pb  is  presented  for  each 
modulation  scheme. 

In  order  to  compare  various  modulation  schemes  of  an  uncoded  2x2  MIMO 
system,  the  average  energy  per  bit  is  calculated  for  each  modulation  scheme  with  regard 
to  the  total  energy  received  during  one  symbol  period. 
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(3.5) 


Let  z  denote  the  received  signal  in  the  absence  of  noise,  given  by 


z 


m 


Therefore,  the  average  total  energy  received  over  one  symbol  period  is 


Et 


m= 1 


(3-6) 


where  fs[»]  is  the  expectation  operator  and  |zm|2  =  zmz* .  Expanding  Equation  (3.6)  we 
have 


N. 


Er  =  ZE 


m=  1 


nt  (  nt  \ 

T.  hm„S„  ^hmkSk 

\k= 1  ) 


n= 1 


nr  ay  nt 

m= 1  n= 1  k= 1 


Due  to  statistical  independence,  Equation  (3.7)  can  be  simplified 

_  Nr  Nt  AY 

ET  hmnSn{hmkSk ) 


m= 1  «=1  £=1 
V*  AY  r 

=11* 

m=l  n= 1 


(3.7) 


(3.8) 


Recall  the  previously  stated  assumptions, 


E 

E 


J  TV  N 


Nr 


(3.9) 


Therefore  the  total  energy  received  during  one  symbol  period  for  a  2x2  MIMO  system 


is 


Mr  IV  T 


Er=ZZE  K 


m= 1  w=l 
Nr  NT  (  , 

yy  — 


f  ,  A 
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(3.10) 
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1. 


BPSK 


The  bit  error  probability  for  BPSK  in  a  single-input  single-output  (SISO)  AWGN 
channel  is  given  by 


Pb=Q 


f  [2R  A 


xtKj 


where  Eb  denotes  the  energy  per  bit  and  Gaussian  Q -function  is  defined  as 

oo 

Q(x)=\-fT=e~y/2dy 

Alternatively,  the  complementary  error  function  is  defined  as 

2 


erfc  (x)  =  —j=  f  e  '  dt 


and  is  related  to  the  Q -function  by 


£>(x)  =  ^-erfc 


vV2  j 


Thus,  the  BER  of  BPSK  in  terms  of  the  complementary  error  function  is 

1 


P,  =  —  erfc 


^  \e  ^ 

c b 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


which  is  more  conducive  to  programming  in  MATLAB. 


Applying  Equation  (3.2)  to  the  case  of  BPSK  modulation,  where  one  bit  of 
information  is  sent  every  symbol,  the  energy  per  bit  is  given  as 


(3.16) 


Figure  13.  shows  the  performance  of  BPSK  in  AWGN,  Rayleigh,  and  Rician  channels 
with  K  -  1  and  K  =  4.  The  theoretical  result  obtained  via  Equation  (3.15)  is  also  plotted, 
which  closely  follows  the  simulated  AWGN  channel. 
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Figure  13.  Simulated  BER  performance  of  a  2x2  MIMO  system 

with  BPSK  modulation. 


2.  QPSK 


The  bit  error  probability  for  QPSK  in  a  SISO  AWGN  channel  is  equal  to  that  of 
BPSK  and  is  given  by 


Pb=Q 


2  Eu 


(3.17) 


where  Eb  denotes  the  energy  per  bit.  Applying  Equation  (3.2)  for  QPSK  modulation  in  a 


2x2  MIMO  system,  where  two  bits  of  infonnation  are  sent  every  symbol,  the  energy  per 
bit  is  given  as 
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1 


(3.18) 


Figure  14.  shows  the  performance  of  QPSK  in  AWGN,  Rayleigh,  and  Rician 
channels  with  K  =  1  and  K  =  4.  The  theoretical  result  obtained  via  Equation  (3.17)  is 
also  plotted.  As  expected,  the  results  for  QPSK  match  those  for  BPSK.  Therefore, 


BPSK  will  not  be  considered  in  the  remaining  analysis. 


Figure  14.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  QPSK 

modulation. 
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3. 


16QAM 


Consider  a  square  16QAM  constellation  depicted  in  Figure  11.  (c).  The  symbol 
error  probability  for  MQAM  in  an  AWGN  channel  when  each  symbol  selected  with 
equal  probability  is  given  by  [2] 


P  =  4 


1- 


1 


Vm  . 


0 


3  log  2M 


r  Eu  3 


M- 1 


(3.19) 


v  '  "  "  "  j 

where  Eb  denotes  the  average  energy  per  bit.  For  the  case  of  16QAM,  Equation  (3.19) 


becomes 


PS  =  1Q 


KNo  j 


(3.20) 


By  utilizing  gray  encoding,  P  can  be  approximated  in  terms  of  Pb  by  [6] 


Pb 


log  2M 


(3.21) 


Applying  Equation  (3.14)  and  Equation  (3.21),  the  BER  of  16QAM  in  terms  of  the 
complementary  error  function, 


-erfc 
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1 

1  Eb 

V 

UJ 

1  N0 

V  ’ 

J 

(3.22) 


The  signal  space  of  a  square  16QAM  constellation  is  depicted  in  Figure  11.  (c)  where  the 
average  symbol  energy  is  given  by 

_  4(2) +  8(10) +  4(1 8) 

16  (3-23) 

=  10 

when  the  minimum  distance  between  adjacent  vectors  is  2  [2],  The  relationship  between 
Eb  and  Es  is 

Eh=^f-  (3.24) 

k 

where  k  denotes  the  bits  per  symbol,  k  =  log2  M  .  Thus,  the  average  energy  per  bit  in 
16QAM  is  found  to  be 


37 


(3.25) 


For  the  case  of  a  2x2  MIMO  system,  we  have 


Eb  = 


5 

2  Nt 


(3.26) 


Figure  15.  shows  the  performance  of  16QAM  in  AWGN,  Rayleigh,  and  Rician 
channels  with  K  =  1  and  K  =  4.  The  theoretical  result  obtained  via  Equation  (3.22)  is 


also  plotted  for  comparison. 


Figure  15.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  16QAM 

modulation. 
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4.  64QAM 


For  the  square  64QAM  constellation  depicted  in  Figure  11.  (d),  the  symbol  error 
probability  is  given  by  [2] 


P  =-Q 
s  2^ 


( 


f  17  \ 


V  N°J 


(3.27) 


v  '  v  '  J 

Applying  Equation  (3.14)  and  Equation  (3.27),  the  BER  of  64QAM  in  terms  of  the 
complementary  error  function  is 


P.  =  — erfc 
4  24 


M3 


v-v 


(3.28) 


and  the  average  symbol  energy  given  by 

4(2)  +  8(10)  +  4(18)  +  8(26)  +  8(34)  +  12(50)  +  8(58)  +  8(74)  +  4(98) 


E  =- 


64 


(3.29) 


=  42 


when  the  minimum  distance  between  adjacent  vectors  is  2  [2].  Using  Equation  (3.24), 
the  average  energy  per  bit  of  64QAM  modulation  is  found  to  be 

Eh  =7  (3.30) 

avg 

For  the  case  of  a  2x2  MIMO  system,  we  have 


Nt 


(3.31) 


Figure  16.  shows  the  perfonnance  of  64QAM  in  AWGN,  Rayleigh,  and  Rician  channels. 
For  comparison,  the  theoretical  result  obtained  via  Equation  (3.28)  is  also  plotted. 
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Uncoded  64QAM 


Figure  16.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  64QAM 

modulation. 

C.  SPACE-TIME  CODED  NARROWBAND 

Building  upon  the  prior  section,  we  evaluate  the  performance  of  space  time 
coding  with  a  MIMO  system.  As  before,  QPSK,  16QAM  and  64QAM  modulations  are 
simulated  and  presented  later  in  this  section.  The  coding  implemented  follows  the 
Alamouti  scheme,  introduced  in  the  previous  chapter  and  the  following  analysis 
completes  that  introduction  by  developing  the  maximum  likelihood  decision  statistic  for 
equal  and  non-equal  energy  constellations. 

Figure  17.  outlines  the  system  model  and  the  output  of  the  combiner,  presented  in 
the  previous  chapter  and  reproduced  here  for  convenience. 
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(3.32) 


+M li  (0  +  Mi  0  +  T)*  +  h2*r \2  (t)  +  M2  ( t  +  T )* 

S  =  (|*n  f  +  1*12  f  +  1*21  f  +  1*22  T  )  ^2 

-Ml  (t  +  T)*  +  *,2*11!  (/)  -  M2  (t  +  T)*  +  h2l\  (t  +  T) 


ih('+r> 


Figure  17.  Alamouti  space-time  coding  scheme  (After  [14]). 


The  Alamouti  scheme  uses  the  following  maximum  likelihood  decision  rule  to 
select  si  iff  [14] 

d2  (y\A\si)  +  d2  [ylAi^d2  ( y\,husk)  +  d 2  ( vf,M,)  (3.33) 

for  all  /'  +  k  .  The  squared  Euclidean  distance  between  signals  x  and  y  is  given  by  [14] 

d2(x,y)  =  (x-y)(x* -y).  (3.34) 

Therefore,  applying  Equations  (3.34)  and  (3.32)  to  the  decision  rule  gives 

( i*i  i  r +i*n  r +1*21  r +1*22  r  -^kt +^2  pu) 

*  (1*1 1 12  +  1*12  I2  +  1*21 12  +  1*22  I2  -  l)  W2  +  d2  (R  +  5,  ) 

(3.35) 

(i*i  i  r + 1*12  r + 1*21  r + 1*22  r  -  i)  k  12 + ^2  (r + ^ ) 

*  (i*i ,  i2 + 1*,2 12 + 1*21 12 + r  - 1)  k  r + d2  (p2 + ) 
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Expanding  Equation  (3.35),  the  maximum  likelihood  statistic  becomes 


mm 


mm 
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(|*„f+|  *, 

(V'nf  +  V’: 
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2  7  2 

12  +  fhi  + 


i/,22r-i)i,t|2-(,*p2+,p2 ) 


where  the 


and 


tenns  of  the  minimum  distance  calculation  are  constant  and  are 


neglected  since  they  have  not  impact  on  the  decision.  For  signals  that  have  equal  energy, 
as  in  the  case  of  BPSK  and  QPSK,  the  magnitudes  are  constant  and  the  statistic  is 


mm 


min 


(3.37) 


-(*,•?, +s,p; 


Figure  18.  shows  the  perfonnance  of  QPSK,  16QAM  and  64QAM  respectively  with 
space  time  coding  in  Rayleigh  and  Rician  channels. 
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Figure  18.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  Alamouti  space- 

time  coding  and  QPSK  modulation. 
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Figure  19.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  Alamouti  space- 

time  coding  and  16QAM  modulation. 
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Space  Time  Coded  64QAM 


Figure  20.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  Alamouti  space- 

time  coding  and  64QAM  modulation. 

D.  SPACE-TIME  CODED  ORTHOGONAL  FREQUENCY  DIVISION 
MULTIPLEXING 

We  now  apply  space-time  coding  to  MIMO-OFDM.  As  a  transmission  strategy, 
MIMO  and  OFDM  are  integral  to  4th  generation  communications  standards.  Thus,  the  bit 
error  performance  of  such  systems  is  of  great  interest.  This  discussion  will  focus  on  the 
implementation  of  STC,  relative  to  the  OFDM  system  presented  in  the  previous  chapter 
[1]. 

Consider  the  basic  space-time  coded  MIMO-OFDM  system  as  shown  in  Figure 
21.  The  coding  scheme  is  the  same  as  detailed  in  the  previous  section;  however,  the 
scheme  is  implemented  in  the  frequency  domain  on  consecutive  tones  or  subcarriers. 
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Analogous  to  the  time  domain  implementation,  in  the  frequency  domain  it  is  assumed 
that  the  channel  remains  constant  over  consecutive  tones  [10]. 


frequency  domain 


Figure  2 1 .  Basic  implementation  of  STC  with  MIMO-OFDM  (After  [5]). 

The  simulation  was  performed  with  an  FFT  size  of  512  points  and  a  cyclic  prefix 
of  1/4.  That  is,  the  last  128  points  from  the  output  of  the  IDFT  operation  were 
reproduced  and  appended  to  the  beginning  of  the  symbol  in  the  form  of  a  cyclic  prefix. 
The  OFDM  symbols  are  then  sent  to  the  transmitter  and  transmitted  through  either  a 
Rayleigh  or  Rician  fading  channel.  At  the  receiver,  the  cyclic  prefix  is  removed  and  the 
DFT  operation  is  performed.  The  symbols  are  then  combined  and  estimated  in 
accordance  with  the  Alamouti  scheme  where  the  decision  criterion  is  analogous  to  the 
narrowband  case  presented  in  the  previous  section. 

The  simulation  evaluated  over  two  million  symbols  with  Eb  / Na  increasing  from  0 

to  30dB  in  increments  of  two  dB.  The  results  of  QPSK,  16QAM  and  64QAM  can  be 
seen  in  Figure  22.  Figure  23.  and  Figure  24.  respectively. 
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Figure  22.  Simulated  BER  performance  of  a  2x2  MIMO-OFDM  system  with  Alamouti 

space-time  coding  and  QPSK  modulation. 
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Figure  23.  Simulated  BER  performance  of  a  2x2  MIMO-OFDM  system  with  Alamouti 

space-time  coding  and  16  QAM  modulation 
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Space  Time  Coded  OFDM,  64QAM 


Figure  24.  Simulated  BER  performance  of  a  MIMO-OFDM  system  with  Alamouti 

space-time  coding  and  64QAM  modulation. 

E.  NARROWBAND  WITH  CHANNEL  STATE  INFORMATION  AT  THE 
TRANSMITTER 

This  section  investigates  the  perfonnance  of  a  2x2  MIMO  system  when  the 
channel  state  is  available  to  the  transmitter,  often  referred  to  as  a  closed-loop  system.  We 
assume  that  the  antennas  are  fixed  and  that  the  feedback  channel  provides  perfect 
information  on  the  transmit  channel.  Figure  25.  depicts  the  basic  system.  Knowledge  of 
the  channel  gains  allows  the  transmitter  to  use  various  techniques  in  order  to  increase 
capacity  or  reliability.  The  following  analysis  focuses  on  increasing  reliability  through 
diversity  techniques.  One  such  method,  which  requires  channel  state  information  at  the 
transmitter,  is  dominate  eigenmode  transmission. 
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Dominate  eigenmode  transmission  is  also  known  as  linear  diversity  precoding, 
and  is  analogous  to  maximal  ratio  Combining  (MRC)  in  the  case  of  a  MISO  system. 
Additionally,  dominate  eigenmode  transmission  lends  a  comparison  with  the  Alamouti 
scheme  of  space-time  coding,  which  does  not  require  CSIT.  Both  techniques  provide  full 
diversity  order  of  four,  the  maximum  in  the  case  of  a  2x2  system.  However,  in  the  case 
of  dominate  eigenmode  transmission,  the  channel  state  information  allows  the  transmitter 
to  optimally  precode  the  signal,  which  provides  an  array  gain  at  the  receiver  not  found  in 
the  Alamouti  scheme. 


feedback 


Figure  25.  Closed  loop  model. 

Given  a  2x2  MIMO  system,  let  the  same  symbol  be  transmitted  on  both 
antennas  during  each  time  period,  weighted  by  a  NT  x  1  vector  w .  Thus,  the  receiver 
yields 

y  =  Hws  +  r|  (3.38) 

where  H  represents  the  NR  x  Nr  channel  matrix,  p  is  the  N R  x  1  noise  vector,  s 
represents  the  transmitted  symbol  and  y  is  the  NR  x  1  received  signal  vector.  At  the 
receiver,  the  signals  are  combined  and  weighted  by  a  NR  x  1  vector  g ,  given  by  [16] 

z  =  gHy  (3.39) 

resulting  in 
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z  =  gH  (Hws  +  ri) 


(3.40) 


In  order  to  find  the  optimal  weighting  of  the  transmitted  signals,  the  transmitter 
utilizes  the  singular  values  of  the  channel.  Recall  from  Chapter  II  that  SVD  operation  is 
given  as 

H  =  USV^  (3.41) 

where  U  and  V  are  Nr  x  Nr  and  NR  x  N R  unitary  matrices,  and  X  is  a  NR  x  NT  diagonal 
matrix  of  the  singular  values  of  the  channel.  It  is  shown  in  [4]  that  the  receiver  SNR  is 
maximized  when  both  the  postcode  vector  g ,  and  the  precode  vector  w  correspond  to  the 
right  and  left  singular  vectors  of  the  maximum  singular  value  of  H ,  crmax .  Therefore,  to 
maximize  the  receiver  SNR,  Equation  (3.40)  can  be  represented  as 

y  =  (7mxs  +  rl  (3-42) 


where  crmax  is  a  scalar  that  corresponds  to  the  largest  value  in  X  which  is  the  maximum 
singular  value  of  H .  Given  that  the  crmax  singular  value  of  H  corresponds  to  the  left  and 
right  singular  vectors  wmax  and  vmax  respectively,  we  can  express  Equation  (3.42)  in  tenns 
of  the  SVD  operation  by  [17] 

y  =  R,ax"HE™xS  +  ri  (3.43) 


The  corresponding  system  model  is  depicted  in  Figure  26.  Thus,  with  channel 
state  infonnation,  the  transmitter  weighs  the  transmitted  symbols  according  to  the  right 
singular  vector  of  the  channel  H .  At  the  receiver,  the  symbols  are  summed  and  weighed 
by  the  left  singular  vector  where  the  signal  to  noise  ratio  is  given  by  [3] 

r  =  <3-44) 

( 7 


where Es  is  the  symbol  energy,  cr  is  the  noise  power  of  each  element  of  r|  and  y  is 


bounded  by  [4], [9] 


Es  H 
crAG 
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E  H 


cr 


(3.45) 
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and  the  Frobenius  nonn  denoted  as  ||[|  is  given  by  [17] 


(3.46) 


Thus,  dominant  eigenmode  transmission  provides  two  types  of  gain;  diversity  and  array. 
Since  the  same  symbol  is  sent  from  each  antenna,  the  diversity  gain  is  the  result  of  having 
multiple  independent  copies  of  the  symbol  for  decision.  At  the  receiver  the  optimally 
weighted  and  independent  signals  are  coherently  combined,  providing  the  array  gain. 
The  array  gain  realized  depends  on  the  fading  characteristics  of  the  channel  and  is  given 
by  £[<7, ;„]  [3][5]. 
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Figure  26.  Dominant  eigenmode  transmission  in  a  2x2  MIMO  system. 


Figure  27.  provides  the  performance  of  QPSK  utilizing  a  dominant  eigenmode 
transmission  in  Rayleigh  and  Rician  channels. 
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Figure  27.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  QPSK  modulation 
using  dominant  eigenmode  transmission  when  the  transmitter  has  channel  state 

information. 


F.  UNKNOWN  RECEIVER 

Throughout  this  chapter,  each  section  presented  analysis  and  performance  results 
of  various  techniques  used  in  MIMO  communications.  The  implementation  of  such 
schemes  is  often  dependant  on  the  environment  in  which  the  communication  system  is 
operating.  Additionally,  the  many  benefits  of  MIMO  communications  rely  on  the 
successful  pairing  of  the  transmitter  and  receiver,  where  the  result  of  this  pairing  is  often 
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the  optimization  of  the  transmission  link.  This  is  a  preferred,  indeed  desired,  outcome  for 
system  designers.  However,  this  thesis  seeks  to  investigate  how  a  receiver  performs  in 
the  absence  of  this  optimization. 

Consider  a  MIMO  system  containing  transmitter  A  and  receiver  B  as  shown  in 
Figure  28.  The  system  is  located  on  fixed  infrastructure  and  has  been  designed  to 
optimize  the  link  between  A  and  B,  a  NLOS  situation.  The  transmitter  has  channel  state 
information,  through  a  feedback  loop  with  receiver  B.  With  CSIT,  the  transmitter  utilizes 
dominant  eigenmode  transmission  in  order  to  improve  the  signal  to  noise  ratio  at  receiver 
B,  providing  a  robust  communications  link. 


Figure  28.  A  fixed  MIMO  system  optimized  for  channel  AB. 

We  now  introduce  a  second  receiver,  C,  that  is  disadvantaged  in  the  sense  that  it 
does  not  enjoy  the  same  optimization  with  transmitter  A  that  receiver  B  has.  Receiver  C 
can  be  mobile  or  fixed.  Channel  AC  can  be  either  Rayleigh  faded  channel  (NLOS)  or  a 
Rician  fading  channel  (LOS).  If  the  receiver  is  mobile,  the  velocity  is  assumed  to  be  such 
that  the  assumptions  given  at  the  beginning  of  this  chapter  still  hold.  Additionally,  the 
receiver  is  not  limited  by  the  number  of  antennas  available.  That  is,  the  receiver  may 
employ  additional  antennas,  greater  than  two,  in  order  to  improve  reception.  Figure  29. 
provides  the  system  layout,  where  the  dashed  line  represents  the  link  that  is  not 
optimized. 
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Figure  29.  Disadvantaged  receiver  in  a  fixed  MIMO  system  optimized  for  channel  AB. 


Recall  from  the  previous  section  that  a  transmitter  with  channel  state  information 
may  precode  the  transmitted  symbols  in  order  to  increase  the  signal  to  noise  ratio  at  the 
receiver.  This  result  was  shown  in  Equation  (3.42)  and  is  reproduced  here  for 
convenience. 

y=crmaxs+'n_  (3-47) 

Expressing  Equation  (3.47)  in  terms  of  the  SVD  operation  we  get  an  equivalent  system 
represented  by  [17] 

y  =  (CxHE™xs  +  fi)^  (3-48) 

where  (□)  denotes  the  characteristics  of  the  link  from  transmitter  A  to  receiver  B.  We 
assume  that  receiver  C  has  knowledge  of  the  precoding  vector  or  the  right  singular  value 
of  the  channel  (H)  .  With  this  knowledge,  the  receiver  is  able  to  compensate  for  the 


optimization  induced  by  the  transmitter  for  channel  AB.  This  case  is  simulated  with 
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QPSK  modulation  over  four  million  symbols  with  Eb/Na  increasing  by  increments  of 
two  dB.  Three  situations  are  presented:  The  first  involves  a  Rayleigh  channel,  the  second 
involves  a  Rician  channel  with  rice  factor  of  one  [K  =  l),  and  the  final  situation  simulates 

a  Rician  channel  with  (if  =  4) .  Additionally,  all  three  scenarios  are  evaluated  with 

receiver  C  having  2,  4,  and  8  antennas.  The  results  are  plotted  in  Figure  30.  Figure  31. 
and  Figure  32.  „  respectively. 


Figure  30.  Simulated  BER  performance  of  the  disadvantaged  receiver  in  a  Rayleigh 

channel,  for  2x2,  2x4,  and  2x8  MIMO  systems. 
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Figure  3 1 .  Simulated  BER  performance  of  the  disadvantaged  receiver  in  a  Rician  channel 

(K=l),  for  2x2,  2x4  and  2x8  MIMO  systems. 
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Figure  32.  Simulated  BER  performance  of  the  disadvantaged  receiver  in  a  Rician  channel 

(K=4),  for  2x2,  2x4,  and  2x8  MIMO  systems. 


G.  SUMMARY 

Chapter  III  provided  bit  error  rate  analysis  and  simulation  results  of  the  concepts 
introduced  in  Chapter  II.  Simulation  assumptions  were  presented  and  used  throughout 
the  chapter  in  an  effort  to  provide  for  results  that  will  be  compared  in  the  next  chapter. 
The  basic  uncoded  narrowband  MIMO  system  was  analyzed  and  simulated  for  numerous 
channel  conditions  as  the  basis  for  follow  on  analysis.  Next,  space  time  coding  following 
the  Alamouti  scheme  was  refined  and  simulated.  Then  orthogonal  frequency  division 
multiplexing  with  space  time  coding  was  simulated.  Channel  state  information  was  then 
incorporated  into  the  narrowband  system  and  simulated.  Finally,  the  chapter  concludes 
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with  a  situation  in  which  a  disadvantaged  receiver  was  introduced  to  an  existing  system 
previously  optimized  for  communication.  Various  results  from  these  simulations  are 
presented  and  compared  in  the  next  chapter. 


59 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


60 


IV.  COMPARISON  AND  RESULTS 


This  chapter  compares  the  performance  of  the  individual  simulations  presented 
the  in  last  chapter.  The  comparison  begins  by  investigating  the  performance  of  an 
uncoded  system  in  Rayleigh  and  Rician  channels.  A  space-time  coded  narrowband 
system  is  compared  to  an  uncoded  system.  Then,  narrowband  performance  is  compared 
to  that  of  OFDM.  Next,  the  narrowband  system  is  analyzed  with  and  without  channel 
state  information.  Finally,  the  performance  of  the  primary  receiver  is  compared  to  that  of 
the  disadvantaged  receiver. 

A.  RAYLEIGH  AND  RICIAN  CHANNEL  MODEL  COMPARISON 

The  discussion  begins  with  an  evaluation  of  the  performance  of  an  uncoded  2x2 
MIMO  system  in  Rayleigh  and  Rician  fading  channels.  Figure  14.  in  Chapter  III,  section 

B,  shows  the  perfonnance  of  an  uncoded  MIMO  system  over  Rayleigh  and  Rician  fading 
channels.  For  this  comparison  QPSK  was  used,  however,  the  results  are  similar  for  the 
other  modulation  schemes. 

Recall  from  Chapter  II,  that  the  K  -  factor  is  a  ratio  of  the  power  in  the  specular 
component  (LOS)  to  that  of  the  diffuse  component  (NLOS).  Clearly,  as  the  value  of  K 
increases,  the  specular  component  dominates  and  the  channel  approaches  that  of  a  non¬ 
fading  channel  or  AWGN  channel.  This  is  confirmed  by  the  increased  perfonnance  of 
the  system  in  a  Rician  channel  to  that  of  the  Rayleigh  channel.  For  a  BER  rate  of  1CT4 ,  a 
LOS  component  of  equal  power  to  that  of  the  NLOS,  as  seen  by  the  receiver,  results  in  an 
improvement  of  4.4  dB.  Table  3.  highlights  the  Eh/Na  required  to  achieve  a  particular 
BER. 
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Table  3.  Eb  / N o  required  to  achieve  various  BER  for  different  channel  models,  2x2  MIMO 

system  with  QPSK  modulation. 


BER 

Rkun(K=41 

BkunO^l) 

Ravish 

Delta 

103 

03  dB 

11.2  dB 

143  dB 

33  dB 

10* 

10.6  dB 

153  dB 

193  dB 

4.4  dB 

lO5 

12.6  dB 

193  dB 

24.7  dB 

5.2  dB 

10* 

14.8  dB 

23.2  dB 

30  dB 

6.8  dB 

Obviously,  when  the  receiver  has  LOS  to  the  transmitter,  higher  order  modulation 
schemes  may  be  employed  and  still  achieve  the  same  (or  better)  BER.  Figure  33. 
compares  QPSK,  16QAM  and  64QAM  with  increasing  power  in  the  specular  component. 
It  is  evident  that  beyond  a  BER  of  10  4  the  higher  modulation  schemes  outperform  QPSK 
due  to  the  specular  component.  This  implies  that  the  specular  component  improves  BER 
more  than  the  higher  order  modulation  degrades  the  BER. 
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Figure  33.  Simulated  BER  performance  of  various  modulation  schemes  with  different 

channel  characteristics,  2x2  MIMO  system. 

B.  UNCODED  AND  SPACE-TIME  CODED  COMPARISON 

We  now  investigate  the  performance  of  an  uncoded  system  to  that  of  a  space-time 
coded  system.  The  Alamouti  space-time  code  is  designed  to  extract  diversity  when  a 
minimum  of  two  transmitting  antennas  are  present,  therefore  we  would  expect  to  see  a 
diversity  gain  when  compared  to  the  uncoded  case.  Figure  34.  plots  the  BER 
performance  of  coded  and  uncoded  QPSK  in  a  Rayleigh  channel. 

The  perfonnance  improvement  seen  the  in  plot  is  due  to  the  diversity  gain 
provided  by  the  Alamouti  scheme.  Additionally,  the  orthogonality  of  the  coding 
sequence  transforms  the  decision  at  the  receiver  from  a  vector  operation  into  a  scalar 
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operation,  reducing  the  receiver  complexity.  Comparing  the  required  Eh/Na  to  achieve 

a  BER  of  1(T4 ,  we  see  that  coding  improves  performance  by  9.7  dB.  Table  4.  provides 
the  required  Eb/No  to  achieve  a  particular  BER. 


Figure  34.  Simulated  BER  performance  of  uncoded  and  coded  2x2  MIMO  system  with 

QPSK  modulation  in  a  Rayleigh  channel. 
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Table  4.  Eb/NB  required  to  achieve  a  certain  BER  in  a  Rayleigh  channel,  2x2  MIMO 
system  using  QPSK  modulation,  with  and  without  Alamouti  space-time  coding. 


men 

HazUgb 

t'nMwtotZmtn 

10* 

14.5  dB 

7.1  dB 

7.4  dB 

10* 

19.8  dB 

lOJdB 

9.7  dB 

10* 

24.7  dB 

12*8  dB 

11.9  dB 

10* 

30  dB 

16  dB 

14  dB 

Next,  we  compare  a  coded  system  to  an  uncoded  system  characterized  by  a  strong 
LOS  component  for  all  three  modulation  schemes,  shown  in  Figure  35.  It  is  interesting  to 
see  from  the  plot  that  an  uncoded  system  over  a  Rician  channel  has  a  marginal 
improvement  in  performance  past  12  dB  (QPSK)  to  that  of  a  space-time  coded  system. 
However,  it  is  clear  that  a  coded  system  outperforms  an  uncoded  system  even  in  an 
environment  where  the  antenna  has  a  dominate  line  of  sight. 

This  demonstrates  the  advantage  of  antenna  placement  with  respect  to  system 
performance,  or  if  this  is  not  possible  the  ability  to  apply  the  Alamouti  scheme  in  order  to 
improve  performance.  MIMO  systems  allow  the  system  designer  implement  multiple 
techniques  in  order  to  establish  robust  communications.  This  simulation  provides  a 
comparison  of  the  perfonnance  of  these  different  techniques.  As  the  simulation  shows, 
Rayleigh  fading  degrades  BER  as  compared  to  Rician  fading  with  a  strong  LOS 
component  (K  =  4)  roughly  the  same  as  uncoded  compared  to  Alamouti  coded. 
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Figure  35.  Simulated  BER  performance  of  a  2x2  MIMO  system,  various  modulation 
schemes,  uncoded  and  using  Alamouti  space-time  coding,  with  different  channel 

characteristics. 

C.  CODED  NARROWBAND  AND  WIDEBAND  COMPARISON 

For  this  thesis,  it  was  assumed  that  the  channel  provided  slow  frequency  non- 
selective  fading  conditions.  Thus,  the  performance  of  a  coded  wideband  system  will 
coincide  with  that  of  a  narrowband  system.  Figure  36.  confirms  this  result  with  a  plot  of 
BER  for  64QAM  with  Alamouti  space-time  coded  wideband  and  narrowband  2x2 
MIMO  system. 
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Figure  36.  Simulated  BER  performance  of  a  2x2  MIMO  system  using  Alamouti  space- 
time  coding,  wideband  and  narrowband  in  various  channels  with  64QAM 

modulation. 

D.  TRANSMITTER  WITH  AND  WITHOUT  CSIT  COMPARISON 

In  this  section,  the  performance  of  a  system  when  the  transmitter  has  CSIT  will  be 
compared  to  a  system  without.  With  a  closed  loop  system,  the  information  available  to 
the  transmitter  can  be  applied  to  improve  performance  in  a  number  of  areas,  depending 
on  the  design  of  the  system.  In  this  thesis  we  focus  on  the  ability  of  the  transmitter  to 
increase  the  received  SNR  through  the  use  of  dominant  eigenmode  transmission,  where 
the  result  of  this  application  is  an  array  gain  in  addition  to  the  diversity  gain.  Therefore, 
we  expect  that  the  performance  increase  over  the  uncoded  case  will  be  considerable. 
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Figure  37.  plots  the  error  rate  performance  of  a  2x2  MIMO  system  with  and 
without  CSIT  where  both  Rayleigh  and  Rician  fading  are  considered.  From  the 
illustration  it  is  clear  that  there  is  a  significant  improvement  when  the  transmitter  has 
knowledge  of  the  channel.  For  a  BER  rate  of  10  4 ,  the  system  with  CSIT  results  in  an 
improvement  of  roughly  12  dB.  Table  5.  highlights  the  Eb/Na  required  to  achieve  a 
particular  BER. 


Figure  37.  Simulated  BER  performance  of  a  transmitter  with  and  without  CSIT  for  a  2x2 

MIMO  system  in  various  channels  using  QPSK  modulation. 
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Table  5.  Eb  /  No  required  to  achieve  a  BER  in  Rayleigh  channel  and  Rician  channels  when 

the  system  has  CSIT,  2x2  MIMO  system. 


BER 

Rayleigh 

CSIT 

Ravkteh 

Delta 

Rician 

CSIT 

Rician 

Delta 

Ifr3 

14.5  dB 

5  dB 

9.5  dB 

11.2  dB 

1.4  dB 

9.8  dB 

HE4 

19.8  dB 

7  3  dB 

11.9  dB 

15.4  dB 

3.8  dB 

11.6  dB 

10* 

24.7  dB 

10.7  dB 

14  dB 

19.5  dB 

63  dB 

133  dB 

10-* 

30  dB 

14  dB 

16  dB 

23.2  dB 

9.5  dB 

13.7  dB 

The  performance  trends  displayed  in  Figure  37.  are  as  expected;  the  array  gain 
increases  the  BER  of  the  entire  system  for  all  values  of  Eb/Na  ,  while  the  diversity  gain 

improves  BER  performance  for  large  values  of  Eb/Na  . 

Figure  38.  compares  a  system  employing  the  Alamouti  space-time  code  and  no 
CSIT  to  a  system  with  CSIT  and  no  STC  in  Rayleigh  fading.  In  this  figure,  it  is  clear  that 
both  systems  offer  the  same  diversity.  However,  the  system  employing  CSIT  shows  a 
consistent  improvement  of  2  dB.  This  is  due  to  an  array  gain  resulting  from  the  coherent 
combining  of  the  precoded  signals  at  the  receiver.  Thus,  for  the  same  diversity  order  we 
see  that  knowledge  of  the  channel  at  the  transmitter  improves  performance  by  2  dB.  Of 
course,  this  improvement  in  perfonnance  is  at  the  cost  of  increased  overhead  due  to  the 
requirements  of  the  feedback  channel. 
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Figure  38.  Simulated  BER  performance  of  a  2x2  MIMO  system  with  CSIT  to  that  of  a 
system  with  Alamouti  space-time  coding  both  in  a  Rayleigh  channel  using  QPSK 

modulation. 


E.  DISADVANTAGED  RECEIVER  COMPARISON 

The  results  in  this  thesis  illustrate  that  MIMO  communications  allow  the  system 
designer  the  ability  to  produce  a  system  that  specifically  addresses  the  environment  in 
which  it  operates.  That  is  one  of  the  benefits  of  MIMO  communications.  When  a  MIMO 
system  is  optimized  using  the  techniques  presented  in  this  thesis,  such  as  employing 
channel  state  information  to  overcome  channel  characteristics  specific  to  that  radio  link, 
introducing  a  second  receiver  to  the  system  could  pose  a  serious  problem,  because  that 
optimization  is  specific  for  the  first  receiver’s  channel  gains. 
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In  this  final  section,  we  investigate  the  performance  of  a  receiver  that  is 
disadvantaged  with  respect  to  the  intended  receiver.  That  is,  the  disadvantaged  receiver 
does  not  enjoy  the  same  optimization  of  the  transmission,  as  does  the  intended  receiver. 
The  goal  of  this  comparison  is  to  integrate  the  topics  discussed  previously  in  order  to 
produce  performance  at  the  disadvantaged  receiver  that  is  equal  if  not  better  than  that  of 
the  intended  receiver. 


Recall,  the  system  in  question,  presented  in  the  last  chapter,  shown  in  Figure  39. 


Figure  39.  Disadvantaged  receiver  in  a  fixed  MIMO  system  optimized  for  channel  AB. 


Transmitter  A  has  channel  state  information  of  channel  AB  and  uses  dominate 
eigenmode  transmission  in  order  to  improve  the  SNR  at  receiver  B,  the  intended  receiver. 
Consider  receiver  C,  which  is  introduced  to  the  system  and  does  not  enjoy  the  same 
optimization  as  the  intended  receiver,  B. 
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1.  Disadvantaged  Receiver  in  a  Rayleigh  Channel 


First,  we  compare  the  performance  of  both  receivers  in  a  Rayleigh  channel,  shown 
in  Figure  40.  The  modulation  for  all  systems  in  this  section  is  QPSK. 


Figure  40.  Simulated  BER  performance  of  a  disadvantaged  receiver  in  a  fixed  MIMO 

system,  Rayleigh  channel. 

From  the  simulation  results,  it  is  clear  that  even  when  the  disadvantaged  receiver 
employs  eight  receive  antennas,  the  performance  is  still  worse  than  for  the  intended 
receiver.  At  a  BER  of  1 0  6  with  eight  receive  antennas,  the  difference  is  1.2  dB. 
Flowever,  depending  on  the  application,  this  performance  may  prove  to  be  adequate. 
Table  6.  provides  performance  values  for  specific  values  of  Eb/N0  . 
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Table  6.  Eh  / No  comparison  for  select  values  of  BER,  disadvantaged  receiver  with  various 

antenna  configurations  over  a  Rayleigh  fading  channel. 


BER 

Bar toteh 

rarntnntatwH 

BajMdkNr=9 

DaHi 

10* 

SdB 

UdB 

13  dB 

10* 

7.9  dB 

lOJdB 

19  dB 

10.7  dB 

lUdB 

15  dB 

10* 

14  dB 

154  dB 

UdB 

2.  Disadvantaged  Receiver  in  a  Rician  Channel  (K  =  1) 

Next,  we  look  at  a  situation  in  which  receiver  C  has  a  LOS  path  to  the  transmitter. 
Where  there  is  equal  power  in  the  specular  and  diffuse  paths,  shown  in  Figure  41. 
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Figure  41 .  Simulated  BER  performance  of  a  disadvantaged  receiver  in  a  fixed  MIMO 

system,  Rician  channel  (K  =  1). 


Clearly,  performance  is  improved  for  receiver  C  compared  to  that  of  the  Rayleigh 
case  of  Figure  40.  ,  as  expected.  From  the  graph,  we  see  that  with  two  receive  antennas 
performance  still  does  not  approach  that  of  the  intended  receiver.  However,  when  the 
number  of  antennas  is  doubled,  the  perfonnance  of  receiver  C  is  comparable  at  low 
values  of  Eb/Na  and  provides  superior  performance  at  high  values  of  Eh/Na  .  Table  7. 
highlights  this  improvement. 
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Table  7.  Eh/Na  comparison  for  select  values  of  BER,  disadvantaged  receiver  with  various 

antenna  configurations  over  a  Rician  channel  (K  =  1). 


BER 

Rjrr toteh 

IMaarfwantaoaH 

Blw.Wr=4 

Petti 

10* 

SdB 

33  dB 

1.1MB 

10* 

7.9  dB 

SJdB 

2.1MB 

10.7  dB 

7.5  dB 

33MB 

10* 

14  dB 

9A  dB 

4.6MB 

3.  Disadvantaged  Receiver  in  a  Rician  Channel  (K  =  4) 

Finally,  we  look  at  a  situation  in  which  receiver  C  has  a  more  significant  LOS 
path  to  the  transmitter,  illustrated  in  Figure  42. 
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Figure  42.  BER  performance  of  a  disadvantaged  receiver  in  a  fixed  MIMO  system, 

Rician  channel  (K  =  4). 

We  observe  that  with  two  antennas,  receiver  C  enjoys  similar  performance  to  the 
intended  receiver.  Obviously,  four  antennas  provide  an  even  greater  performance  gain. 
A  transmitter  with  eight  antennas  was  not  evaluated  since  the  result  of  four  antennas  is 
sufficient,  but  it  is  expected  that  BER  performance  increases  dramatically  as  the  number 
of  antennas  increases.  Table  8.  gives  specific  values  for  Eb/Na  . 
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Table  8.  Eh  / No  comparison  for  select  values  of  BER,  disadvantaged  receiver  with  various 

antenna  configurations  over  a  Rician  channel  (K  =  1). 


BCR 

Intended 

Ravleteh 

Dbidvantased 

Rice.  Nr =4 

Delta 

Dbadvantased 

Rice.  Nr  =  2 

Delta 

1<H 

SdB 

2.5  dB 

2.5  dB 

5j6dB 

6.6  dB 

i<r* 

7.9  dB 

4.2  dB 

3.7  IB 

7.8  dB 

0.1  dB 

KT5 

10.7  dB 

5.1  dB 

5.6  dB 

10.3  dB 

0.4  dB 

1<T* 

14  dB 

6dB 

SdB 

12  dB 

6dB 

When  considering  the  performance  of  the  disadvantaged  receiver  as  the  number 
receive  antennas  increased,  recall  from  the  first  section  in  Chapter  III,  that  the  average 
total  energy  per  bit  received  is  equal  to  the  transmitted  energy  per  bit  per  transmit 
antenna  by  convention.  In  other  words,  both  systems  have  the  same  total  average  Eb .  As 

a  result,  increasing  the  number  of  receive  antennas  reduces  the  channel  gains  of  the 
system.  Thus,  our  assumption  diminishes  the  perfonnance  gained  by  increasing  the 
number  of  receive  antennas  at  the  disadvantaged  receiver.  Of  course,  if  the  systems  had 
different  Eb  and  Na ,  those  values  would  be  used  to  make  a  more  accurate  comparison. 
The  gross  path  loss,  or  geometry  with  respect  to  the  transmitter,  may  favor  one  receiver 
over  the  other.  In  this  case,  we  would  see  an  even  greater  improvement  in  performance 
since  the  greater  number  of  receive  antennas  results  in  more  received  power  and  thus  a 
higher^ . 

However,  with  the  assumptions  of  this  thesis,  the  results  of  the  simulations  show 
that  when  the  disadvantaged  receiver  has  a  weak  Eb/Na ,  increasing  the  number  of 
receive  antennas  is  one  method  to  significantly  decrease  the  bit  error  rate. 

F.  SUMMARY 

In  this  chapter,  several  different  scenarios  are  compared  in  terms  of  BER 
performance,  based  on  the  results  of  the  previous  chapter.  The  first  comparison 
investigated  the  performance  of  an  uncoded  MIMO  system  over  Rayleigh  and  Rician 
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channels.  The  resulting  improvement  in  BER  performance  over  a  Rician  channel  was 
presented.  Next,  a  space-time  coded  system  was  compared  to  an  uncoded  system,  where 
the  diversity  gain  provided  by  the  Alamouti  scheme  was  made  apparent.  Then, 
narrowband  performance  was  compared  to  that  of  OFDM.  Subsequently,  a  MIMO 
system  was  analyzed  with  and  without  channel  state  information.  Finally,  the 
performance  of  the  primary  receiver  in  a  MIMO  system  is  quantified  and  compared  to  the 
performance  of  a  disadvantaged  receiver.  The  disadvantaged  receiver  is  then  presented 
with  various  configurations  in  order  to  improve  it  BER  performance.  The  next  chapter 
concludes  this  thesis  and  provides  a  summary  of  results. 


78 


V.  CONCLUSION 


The  goal  of  this  research  was  to  investigate  disadvantaged  receiver  performance 
compared  to  the  intended  receiver  in  a  MIMO  system  with  the  transmission  optimized  for 
the  intended  receiver.  This  was  accomplished  through  a  progressive  development 
process  from  the  simple  uncoded  MIMO  system  to  one  that  employs  channel  state 
information  at  the  transmitter.  Once  the  final  system  was  developed,  a  second, 
disadvantaged  receiver  was  added  and  configured  such  that  the  BER  perfonnance 
matched  or  exceeded  that  of  the  intended  receiver.  Each  system  was  simulated  in 
MATLAB  and  the  BER  performance  was  determined  and  compared. 

A.  SUMMARY 

A  comprehensive  study  of  MIMO  systems  and  technologies  was  presented. 
Specifically,  the  BER  performance  of  the  transmitter  and  receiver  in  several  cases  was 
developed,  to  include  space-time  coding,  OFDM,  and  channel  state  information  at  the 
transmitter.  A  second  receiver  was  then  introduced  into  the  final  system  and  the 
performance  utilizing  the  discussed  technologies  was  investigated. 

The  systems  were  designed  and  simulated  in  MATLAB  for  each  case,  over  both 
Rayleigh  and  Rician  fading  channels.  The  results  obtained  were  presented  in  the  form  of 
BER  performance  curves  and  compared  against  each  other. 

B.  SIGNIFICANT  RESULTS 

The  following  results  are  provided  by  this  thesis. 

The  BER  performance  of  a  MIMO  system  increases  significantly  when  a  LOS 
component  is  present.  This  is  completely  consistent  with  intuition  developed  from  the 
study  of  SISO  systems.  When  equal  power  is  present  in  the  specular  and  diffuse 
components  (K  =  1),  an  improvement  of  6.8  dB  was  observed  for  a  bit  error  rate  of  1CT6 . 
In  addition,  a  system  employing  a  higher  order  modulation  scheme  such  as  64QAM  over 
a  Rician  channel  is  able  to  achieve  equal,  if  not  better,  perfonnance  than  a  system  with 
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QPSK  modulation  over  a  Rayleigh  channel.  This  comparison  was  made  in  order  to 
quantify  the  advantage  of  a  receiver  in  a  Rician  channel  as  compared  to  one  in  a  Rayleigh 
channel.  That  is  important  in  this  work  as  that  might  be  the  case  for  the  disadvantaged 
and  intended  receivers  described  in  Chapter  IV,  section  E. 

In  Chapter  IV,  section  B,  MIMO  systems  utilizing  the  Alamouti  space-time  code 
were  observed  to  improve  the  BER  performance  by  14  dB  at  a  bit  error  rate  of  1CT6  when 
compared  to  uncoded  systems.  Moreover,  in  high  Eb/Na  regions,  a  system 

communicating  over  a  Rician  channel  provided  increased  performance  to  that  of  a  coded 
system  over  a  Rayleigh  channel.  This  comparison  further  demonstrates  the  advantage  of 
antenna  placement  with  respect  to  system  performance.  This  relates  to  the  situation  of  a 
disadvantaged  receiver  since  the  fixed  system  may  be  designed  in  such  a  way  as  to 
exploit  this  result.  In  other  words,  with  MIMO  systems  there  are  multiple  techniques  that 
can  be  implemented  in  order  to  establish  robust  communications.  This  work  primarily 
focused  on  the  situation  in  which  the  transmitter  utilized  channel  state  information  in 
order  to  establish  a  reliable  system.  However,  understanding  of  such  other  techniques 
benefits  the  designer  of  the  disadvantaged  system  as  well. 

When  the  transmitter  utilizes  channel  state  information  to  increase  receiver 
diversity  we  observe  a  significant  performance  increase  over  a  system  that  does  not  have 
CSIT.  For  a  BER  of  1CT6 ,  an  improvement  of  13.7  dB  was  revealed.  Additionally,  it  was 
observed  that  the  system  with  CSIT  provides  an  array  gain  of  2  dB  when  compared  to  the 
system  without  channel  state  information. 

When  a  disadvantaged  receiver  is  added  to  the  system  it  is  seen  that  the  BER 
performance  cannot  match  that  of  the  intended  receiver  when  both  receivers  experience 
Rayleigh  fading,  even  when  the  number  of  reception  antennas  is  increased  to  eight. 
However,  when  the  disadvantaged  receiver  is  located  such  that  it  has  a  line  of  sight 
component  (Rician  fading),  we  observe  that  performance  increases  significantly.  In  the 
case  of  equal  LOS  and  NLOS  components,  a  disadvantaged  receiver  with  four  reception 
antennas  outperforms  the  intended  receiver.  When  the  disadvantaged  receiver  has  a  more 
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significant  LOS  component,  similar  performance  to  that  of  the  intended  receiver  can  be 
found  with  only  two  receive  antennas  at  the  disadvantaged  receiver. 

C.  RECOMMENDATIONS  FOR  FUTURE  WORK 

Four  areas  have  been  identified  for  future  work.  First,  this  work  assumes  that  the 
system  is  located  in  a  fixed  environment.  However,  a  more  realistic  environment  may  be 
investigated  by  relaxing  the  channel  assumptions  to  allow  for  frequency  selective  fading. 

Second,  the  communication  system  can  be  further  refined  with  a  space-time 
coding  scheme  and  OFDM.  Due  to  the  channel  conditions  assumed,  this  thesis 
developed  a  system  that  was  single  carrier  and  uncoded.  If  the  first  recommendation  is 
applied,  then  the  system  can  be  integrated  with  those  techniques  and  the  disadvantaged 
receiver  performance  can  then  be  determined  under  the  new  channel  conditions.  Of 
course,  utilizing  OFDM  with  the  first  recommendation  will  provide  similar  results  to  the 
flat  fading  conditions  considered  in  this  work,  since  each  subcarrier  experiences  flat 
fading  conditions  even  when  transmitted  through  a  frequency  selective  channel. 
However,  a  simulation  can  be  conducted  to  confirm  this  result. 

Third,  in  this  thesis  the  focus  was  on  bit  error  performance.  Future  work  can 
analyze  the  system  from  an  information  theory  standpoint.  That  is,  focus  on  the  capacity 
of  the  various  systems  and  compare  the  capacity  gains  (or  losses).  For  this,  the  channel 
state  infonnation  at  the  transmitter  can  be  utilized  to  increase  capacity  rather  than 
improving  the  receiver  BER  as  was  used  in  this  study. 

Fourth,  the  simulations  performed  in  this  work  discovered  that  as  the  number  of 
receive  antennas  employed  by  the  disadvantaged  receiver  increased,  so  did  the  BER 
performance.  However,  future  work  can  identify  the  point  of  diminishing  returns  with 
respect  to  the  number  of  receivers.  At  what  point  do  the  number  of  antennas  and  the 
corresponding  complexity  outweigh  the  perfonnance  improvement  gained?  Simulations 
involving  an  increased  number  of  antennas  can  easily  determine  such  a  threshold. 
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